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We study proton-(anti)proton collisions at the LHC or Tevatron in the presence of experimental 
restrictions on the hadronic final state and for generic parton momentum fractions. At the scale Q 
of the hard interaction, factorization does not yield standard parton distribution functions (PDFs) 
for the initial state. The measurement restricting the hadronic final state introduces a new scale 
HB <C Q and probes the proton prior to the hard collision. This corresponds to evaluating the PDFs 
at the scale fiB- After the proton is probed, the incoming hard parton is contained in an initial- 
state jet, and the hard collision occurs between partons inside these jets rather than inside protons. 
The proper description of such initial-state jets requires "beam functions". At the scale /is, the 
beam function factorizes into a convolution of calculable Wilson coefficients and PDFs. Below \xb, 
the initial-state evolution is described by the usual PDF evolution which changes x, while above 
Hb it is governed by a different renormalization group evolution that sums double logarithms of 
I-Ib/Q and leaves x fixed. As an example, we prove a factorization theorem for "isolated Drell-Yan", 
pp — > X£ + £~ where X is restricted to have no central jets. We comment on the extension to cases 
where the hadronic final state contains a certain number of isolated central jets. 



I. INTRODUCTION 

Factorization is one of the most basic concepts for un- 
derstanding data from the Tevatron at Fcrmilab and the 
CERN Large Hadron Collider (LHC). For a review of fac- 
torization see Ref. [J. Typically, factorization is viewed 
as the statement that the cross section can be computed 
through a product of probability functions, namely par- 
ton distribution functions (PDFs), describing the prob- 
ability to extract a quark or gluon from the protons 
in the initial state, a perturbative cross section for the 
hard scattering, and a probabilistic description of the fi- 
nal state by a parton shower Monte Carlo or otherwise. 
This factorization is of key importance in the program 
to search for new physics, as new physics is primarily a 
short-distance modification of the hard scattering that 
must be distinguished from the array of QCD interac- 
tions in the initial and final states. Factorization is also 
necessary for controlling QCD effects. For example, the 
momentum distributions of the colliding partons in the 
protons arc nonperturbative, but factorization can imply 
that these are described by universal distributions which 
have been measured in earlier experiments. 

As the primary goal of the experiments at the LHC or 
Tevatron is to probe the physics of the hard interaction, 
measurements often impose restrictions on the hadronic 
final state, requiring a certain number of hard leptons or 
jets in the final state [IHEj]- For example, a typical new 
physics search looking for missing transverse energy may 
also require a minimum number of jets with pt above 
some threshold. To identify the new physics and deter- 
mine the masses of new-physics particles, one has to re- 
construct decay chains with a certain number of jets and 
leptons in the final state. 

Any theoretical prediction for pp or pp collisions, 
whether analytic or via Monte Carlo generators, depends 
on factorization. However, for the majority of processes 
of interest at hadron colliders where one distinguishes 



properties of the hadronic final state, so far no rigor- 
ous field-theoretic derivation of a factorization theorem 
to all orders in perturbation theory exists. The most 
well-known factorization theorem is 

da = ^d ( 7P art ®/ 4 (e a )®/,(6), (1) 

where fi and fj arc the standard PDFs for partons 
hj = {ff: u i d, . . .} carrying momentum fractions £ a 
and £b (which we use as our PDF x- variables), and der^ art 
is the partonic cross section to scatter i and j calcu- 
lated in fixed-order perturbation theory. In Eq. ([1]), the 
hadronic final state is treated as fully inclusive. Hence, 
in the presence of experimental restrictions that make a 
process less inclusive, Eq. dTJ is a priori not applicable. 
At best, an additional resummation of large phase-space 
logarithms must be carried out by a further factorization 
of dcr^ alt , while at worst, additional nonperturbative in- 
formation beyond that contained in the PDFs is required 
or there is no factorization. 

Factorization theorems for threshold resummation in 
hadron-hadron collisions are a well-studied case where 
Eq. (|T|) can be extended to sum large phase-space loga- 
rithms The corresponding formalism however re- 
quires the limit x — ¥ 1, and hence is not directly relevant 
at the LHC, where the cross section for most measure- 
ments is dominated by the region x far from one (l5j . 

Our goal is to study factorization for a situation where 
the hard interaction occurs between partons with generic 
momentum fractions, away from the limit x 1, and 
where the hadronic final state is measured and restricted 
by constraints on certain kinematic variables. These re- 
strictions allow one to probe more details about the final 
state and may be used experimentally to isolate central 
hard jets or leptons or to control backgrounds. 

A typical event at the LHC with three high-p^ jets is 
illustrated in Fig. [1] There are several complications one 
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FIG. 1: A typical event with jet production at the LHC. 



has to face when trying to derive a factorization theo- 
rem in this situation. First, experimentally the number 
and properties of the final-state jets are determined with 
a jet algorithm. Second, to enhance the ratio of signal 
over background, the experimental analyses have to ap- 
ply kinematic selection cuts. Third, in addition to the 
jets produced by the hard interaction, there is soft ra- 
diation everywhere (which is part of what is sometimes 
called the "underlying event" ) . Fourth, a (large) fraction 
of the total energy in the final state is deposited near the 
beam axes at high rapidities. An important component 
of this radiation can contribute to measurements, and 
when it does, it cannot be neglected in the factorization. 
In this paper we focus on the last three items. Methods 
for including jet algorithms in factorization have been 
studied in Refs. [1, 111, 03 

To allow a clean theoretical description, the observ- 
ables used to constrain the events must be chosen care- 
fully such that they are infrared safe and sensitive to 
emissions everywhere in phase space. Observables satis- 
fying these criteria for hadron colliders have been clas- 
sified and studied in Refs. [H, [l9j|, and are referred to 
as global event shapes. (Issues related to non-global ob- 
servables have been discussed for example in Refs. [20l — 
[23j|.) For our analysis we use a very simple example of 
such an observable, constructed as follows. We define two 
hemispheres, a and b, orthogonal to the beam axis and 
two unit lightlike vectors n a and rib along the beam axis 
pointing into each hemisphere. Taking the beam axis 
along the z direction, hemisphere a is defined as z > 
with n£ = (1,0,0,1), and hemisphere b as z < with 
n b = (1,0,0,-1). We now divide the total momentum 
px of the hadronic final state into the contributions from 
particles in each hemisphere, px = Px a + Px b ■ Next, we 
remove the momenta pj of all jets (defined by an appro- 
priate jet algorithm) in each hemisphere. Of the remain- 
ing hemisphere momenta, we measure the components 
.B+ and B b defined by 

BZ=n a -(p Xa -J2Pj), ( 2 ) 

,j£a 



and analogously for B b . Because of the dot product with 
n a or rib, energetic particles near the beam axes only give 
small contributions to B£ or B b . In particular, any con- 
tributions from particles at very large rapidities outside 
the detector reach, including the remnant of unscattered 
partons in the proton, arc negligible. All observed par- 
ticles contribute either to B~£ , B b , or a jet momentum, 
so we are ensured that we cover all of phase space. De- 
manding that B^ b are small restricts the radiation be- 
tween central jets, only allowing highly energetic particles 
either within these jets or inside jets along the beam di- 
rections labeled "Jet a" and "Jet b" in Fig. [TJ Hence, 
measuring and constraining B^ b provides a theoretically 
clean method to control the remaining particles in the 
hadronic final state. This ensures that observables based 
on the large momenta of hard jets or leptons are clean, 
safe from uncontrolled hadronic effects. 

In this paper, we consider the simplest situation where 
the above setup can be realized, allowing us to explore the 
implications of restrictions on the hadronic final state. 
We prove a factorization theorem for Drell-Yan produc- 
tion pp — > X£ + £~ where X is allowed to have hard jets 
close to the beam, but no hard central jets. We call this 
"isolated Drell-Yan" . Our proof of factorization uses the 
soft-collinear effective theory (SCET) plus ad- 

ditional arguments to rule out possible Glauber effects 
based in part on Refs. [28|,[2!|. Although we focus our dis- 
cussion on Drell-Yan, our factorization theorem applies 
to processes pp — > XL, were the lepton pair is replaced 
by other non-strongly interacting particles, such as Higgs 
or Z' decaying non-hadronically. Though our analysis is 
only rigorous for pp — > XL, we also briefly discuss what 
the extended factorization formula may look like for pro- 
cesses with additional identified jets in the final state. 

Our main result is to show that process-independent 
"beam functions", Bi(t,x) with i = {g,u,u,d, . . .}, are 
required to properly describe the initial state. For the 
usual PDFs in Drell-Yan production appearing in Eq. (fTJ) , 
the hadronic final state X is treated fully inclusively, and 
the effects of initial- and final-state soft radiation cancel 
out [l[. With restrictions on X, the effects of soft ra- 
diation can no longer cancel. Generically, by restricting 
X one performs an indirect measurement of the proton 
prior to the hard collision. At this point, the proton is 
resolved into a colliding hard parton inside a cloud of 
collincar and soft radiation. The proper description of 
this initial-state jet is given by a beam function in con- 
junction with an appropriate soft function describing the 
soft radiation in the event. 

One might worry that the collision of partons inside 
initial-state jets rather than partons inside protons could 
drastically change the physical picture. Although the 
changes are not as dramatic, they have important im- 
plications. The beam function can be computed in an 
operator product expansion, giving 

Bi(t, £, MB ) = S(t) f^, n B ) + 0[a a (jjL B )] , (3) 
where hb is an intermediate pcrturbative scale and t is an 
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invariant-mass variable closely related to the off-shellness 
of the colliding parton (and the Mandelstam variable t). 
Thus, the beam functions reduce to standard PDFs at 
leading order. For what we call the gluon beam function, 
this was already found in Ref. [30j , where the same matrix 
element of gluon fields appeared in their computation of 
J/tpX using SCET. 

Equation <j3j> implies that the momentum fractions £ Qi b 
are determined by PDFs evaluated at the scale \xb -C Q, 
which is parametrically smaller than the scale Q of the 
partonic hard interaction. The renormalization group 
evolution (RGE) for the initial state now proceeds in 
two stages. For scales \x < ur, the RGE is given by 
the standard PDF evolution [3ll - l35| . which sums single 
logarithms, mixes the PDFs, and redistributes the mo- 
mentum fractions in the proton to lower x values. For 
scales /i > /ig, the jet- like structure of the initial state 
becomes relevant and its evolution is properly described 
by the RGE of the beam function. In contrast to the 
PDF, the evolution of the beam function is independent 
of x, does not involve any mixing between parton species, 
and sums Sudakov double logarithms. In addition to the 
change in evolution, the transition from PDFs to beam 
functions at the scale \ib also involves explicit u s {[ib) 
corrections as indicated in Eq. ([3]) . These include mixing 
effects, such as a gluon from the proton pair-producing 
a quark that goes on to initiate the hard interaction and 
an antiquark that is radiated into the final state. For our 
observables such fluctuations are not fully accounted for 
by the PDF evolution. These beam effects must be taken 
into account, which can be done by perturbative calcula- 
tions. The standard PDFs are still sufficient to describe 
the nonperturbative information required for the initial 
state. 

One should ask whether the description of the ini- 
tial state by beam functions, as well as their interplay 
with the soft radiation, are properly captured by current 
Monte Carlo event generators used to simulate events at 
the LHC and Tevatron, such as Pythia [HI, H?} and Her- 
wig 38, 39]. In these programs the corresponding effects 
should be described at leading order by the initial-state 
parton shower in conjunction with models for the under- 
lying event [40M43| . The experimental implications and 
reliability of these QCD Monte Carlo models have been 
studied extensively |44l - t46| . We will see that the initial- 
state parton shower is in fact closer to factorization with 
beam functions than to the inclusive factorization for- 
mula in Eq. ([!]). In particular, the physical picture of 
off-shell partons that arises from the factorization with 
beam functions has a nice correspondence with the pic- 
ture adopted for initial-state parton showers a long time 
ago JUS]. There are also differences. Our analysis is 
based solely on QCD soft-collinear factorization, whereas 
the initial-state parton shower is partly based on the pic- 
ture arising from small- a; physics or semihard QCD |49| . 
For the parton distributions our formalism applies in a 
situation that is intermediate between the case of very 
small x, where a resummation of In a; becomes impor- 



tant, and the case x — >■ 1, where threshold resummation 
in ln(l — x) becomes important. Numerically, our results 
apply for the dominant region of x values that arc of in- 
terest at the LHC. Experimentally, measurements of the 
isolated Drell-Yan cross section provide a simple obsev- 
able that can rigorously test the accuracy of the initial- 
state shower in Monte Carlo programs, by contrasting it 
with the analytic results reported here. 

In Sec. [HI we discuss our main results and explain vari- 
ous aspects of the factorization with beam functions. The 
goal of this section is to give a thorough discussion of the 
physical picture behind our results which is nontechni- 
cal and accessible to non-expert readers. In Sec. IIIII we 
elaborate on the field-theoretic definition and properties 
of the beam functions and their relation to the PDFs. 
We quote explicit results for the quark beam function at 
one loop, the derivation of which will be given in a sepa- 
rate publication [5(|. In Sec. IIV[ we derive in detail the 
factorization theorem for isolated pp —> XL using SCET, 
and apply it to the case of Drell-Yan. Readers not inter- 
ested in the technical details can freely skip this section. 
Plots of the isolated Drell-Yan cross section are given in 
Sec. El We conclude in Sec. [VlJ 



II. FACTORIZATION WITH BEAM 
FUNCTIONS 

This section provides an extensive discussion of how 
factorization with beam functions works, including the 
necessary kinematic definitions for the variables that con- 
strain the hadronic final state. In the interest of avoiding 
technical details, we only discuss the physics contained 
in the factorization theorems. Readers interested in the 
field-theoretic definitions for the beam functions are re- 
ferred to Sec. IIIII while those interested in the derivation 
of the factorization theorem in SCET and explicit defini- 
tions for all its ingredients are referred to Sec. IIVI 

In Sec. Ill At we review the factorization theorems for 
inclusive Drell-Yan and threshold Drell-Yan, and then 
explain the factorization theorem for our isolated Drell- 
Yan process. We use a simple setup where measurements 
on the final-state hadrons use hemispheres orthogonal to 
the beam. These observables are generalized in Sec. Ill Bl 
to uniformly account for measurements that sample over 
a wide variety of boosts between the hadronic and par- 
tonic center-of-mass frames. We explain the relation be- 
tween beam functions and parton distribution functions 
in Sec. Ill CI We compare the beam-function renormal- 
ization group evolution to initial-state parton showers in 
Sec. Ill Dl In Sec. Ill El we show how the various pieces 
in the factorization theorem arise from the point of view 
of a fixed-order calculation. In Sec. Ill Ft we compare the 
structure of large logarithms and their resummation for 
the different factorization theorems. This yields an in- 
dependent argument for the necessity of beam functions 
and provides a road map for incorporating beam func- 
tions in other isolated processes. Finally in Sec. Ill Gl we 
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comment on the extension of the factorization with beam 
functions to the case where one has two or more isolated 
jets in the final state. 



A. Drell-Yan Factorization Theorems 

To describe the Drell-Yan process pp — > X£ + £~ or 
pp —> X£ + £~ , we take 



(4) 



where P£ b are the incoming (anti)proton momenta, 

E cm = -J (P a + P ) 2 is the total center-of-mass energy, 
and is the total momentum of the £ + £~ pair. We also 
define 



E 2 



2 P a ■ q 
Xb = Vre~ Y , 



(5) 



where Y is the total rapidity of the leptons with respect 
to the beam axis, and x a and Xb are in one-to-one corre- 
spondence with r and Y. Their kinematic limits are 



< r < 1, 

T < X a < 1 



2\Y\ < -Inr, 

T < Xb < 1 . 



(6) 



The invariant mass of the hadronic final state is bounded 
by 



'x 



= p\ < E 2 m (l - V7) 2 . 



In Drell-Yan 



q 2 > A 



QCD 



(7) 



(8) 



plays the role of the hard interaction scale. In general, 
for factorization to be valid at some leading level of ap- 
proximation with a perturbative computation of the hard 
scattering, the measured observable must be infrared safe 
and insensitive to the details of the hadronic final state. 



For inclusive Drell-Yan, illustrated in Fig. 2(a) 



sums over all hadronic final states X allowed by Eq. ([7|) 
without imposing any cuts. Hence, the measurement is 
insensitive to any details of X because one sums over all 
possibilities. In this situation there is a rigorous deriva- 
tion of the classic factorization theorem |28l k51L k32J 



1 da 



cr dq 2 dY 



E 

i,3 



H 



incl 



(Xa_ 



X b 2 

q ,/•* 



xM£a,M)/;(6,M)[l + o( 



' AqcdN 
Q > 



(9) 



where cto = 4nal m / (3N c E 2 m q 2 ), and the integration lim- 
its are x a < £ a < 1 and Xb < £,b < 1- The sum is 
over partons i,j = {g, u, u,d, . . .}, and fi(£ a ) is the par- 
ton distribution function for finding parton i inside the 
proton with light-cone momentum fraction £ a along the 



proton direction. Note that £ aj {, are partonic variables, 
whereas x a ,b are leptonic, and the two are only equal at 
tree level. The inclusive hard function Hf} can be com- 
puted in fixed-order perturbative QCD as the partonic 
cross section to scatter partons i and j [corresponding to 
dcr,^ alt in Eq. (JTJ)] and is known to two loops [53M57| . 

For threshold Drell-Yan, one imposes strong restric- 
tions to only allow soft hadronic final states with mx <C 
Q, as illustrated in Fig. 2(b) Using Eq. ([7]), this can be 
ensured by forcing (1 — -y/r) z <C r, so that one is close to 
the threshold r — > 1. In this case, there are large dou- 
ble logarithms that are not accounted for by the parton 
distributions. Furthermore, since 



1 > £a,b > x a ,b > t ->• 1 , 



(10) 



a single parton in each proton carries almost all of the 
energy, £ a< b —t 1. The partonic analog of r is the variable 



<t 



£a£,bE 2 m Ca6 



<1, 



(11) 



and r — > 1 implies the partonic threshold limit z — > 1. 
As Eq. Q forces Y — > for r — > 1, it is convenient to 
integrate over Y and consider the r — > 1 limit for da/dc] 2 
The relevant factorization theorem in this limit is [a, 

1 S = E tf> f^r^ /*(&,/*) *) 

a o aq ~f J t,a t,b 

xQS thT [Q 



1 - 



£a6 



)( AQCD i r 



Q 



(12) 



where we view Eq. (fl2|) as a hadronic factorization theo- 
rem in its own right, rather than simply a rcfactorization 
of H'^ cX in Eq. ©. This Drell-Yan threshold limit has 
been studied extensively [l2|, EH [H|-[62j]. Factorization 
theorems of this type are the basis for the resummation of 
large logarithms in near-threshold situations. In contrast 
to Eq. ©, the sum in Eq. (|12[) only includes the domi- 
nant qq terms for various flavors, ij = {uu,uu,dd, . . .}. 
Other combinations are power-suppressed and only ap- 
pear at 0(1 — r) or higher. The threshold hard function 
Hij ~ \CiC* | is given by the square of Wilson coeffi- 
cients in SCET, and can be computed from the timelike 
quark form factor. The threshold Drell-Yan soft function 
Sthr i s defined by a matrix element of Wilson lines and 
contains both perturbative and nonperturbative physics. 
If it is treated purely in perturbation theory at the soft 
scale Q(l — r), there are in principle additional power 
corrections of C[A QCD /Q(1 - r)] in Eq. (TTJ (H. 

Our goal is to describe the isolated Drell-Yan process 
shown in Fig. 2(c) Here, the colliding partons in the hard 
interaction are far from threshold as in the inclusive case, 
but we impose a constraint that does not allow central 
jets. Soft radiation still occurs everywhere, including the 
central region. Away from threshold, the hard interaction 
only carries away a fraction of the total energy in the col- 
lision. The majority of the remaining energy stays near 



5 




(a) Inclusive Drell-Yan production. (b) Drell-Yan near threshold. (c) Isolated Drell-Yan. 




(d) Dijet production near threshold. (e) Isolated dijet production. 



FIG. 2: Different final-state configurations for pp collisions. The top row corresponds to Drell-Yan factorization theorems for 
the (a) inclusive, (b) threshold, and (c) isolated cases. The bottom row shows the corresponding pictures with the lepton pair 
replaced by dijets. 



the beam. The colliding partons emit collinear radiation 
along the beams that can be observed in the final state, 
shown by the green lines labeled "Jet a" and "Jet 6" in 
Fig. 2(c) This radiation cannot be neglected in the fac- 



torization theorem and necessitates the beam functions. 
In the threshold case, these jets are not allowed by the 
limit t — y 1 , which forces all available energy into the 
lcptons and leaves only soft hadronic radiation. 1 In the 
inclusive case there are no restrictions on additional hard 
emissions, in which case initial-state radiation is included 
in the partonic cross section in Hf 



rincl 
ij ' 



Also shown in Fig. 2(c) is the fact that the leptons 
in isolated Drell-Yan need not be back-to-back, though 
they are still back-to-back in the transverse plane [see 
Sec. IIV B| . In this regard, isolated Drell-Yan is in- 



1 Note that the proof of factorization for the partonic cross section 
in the partonic threshold limit z — > 1 is not sufficient to establish 
the factorization of the hadronic cross section, unless one takes 
the limit r — + 1. The hadronic factorization theorem assumes 
that all real radiation is soft with only virtual hard radiation 
in the hard function. The weaker limit z — ¥ 1 still allows the 
incoming partons to emit energetic real radiation that cannot 
be described by the threshold soft function. Only the r — > 1 
limit forces the radiation to be soft. This point is not related to 
whether or not the threshold terms happen to dominate numer- 
ically away from r — > 1 due to the shape of the PDFs or other 



between the threshold case, where the leptons are fully 
back-to-back with Y ~ 0, and the inclusive case, where 
they are unrestricted. 

In Figs. 12(5)1 and HM 

we show analogs of threshold 
Drell-Yan and isolated Drell-Yan where the leptons are 
replaced by final-state jets. We will discuss the extension 
to jets in Sec. Ill Gl below. 

To formulate isolated Drell-Yan we must first discuss 
how to veto hard emissions in the central region. For this 
purpose, it is important to use an observable that covers 
the full phase space. Jet algorithms are good tools to 
identify jets, but not necessarily to veto them. Imagine 
we use a jet algorithm and require that it does not find 
any jets in the central region. Although this procedure 
covers the full phase space, the restrictions it imposes 
on the final state depend in detail on the algorithm and 
its criteria to decide if something is considered a jet or 
not. It is very hard to incorporate such restrictions into 
explicit theoretical calculations, and in particular into a 
rigorous factorization theorem. Even if possible in prin- 
ciple, the resulting beam and soft functions would be 
very complicated objects, and it would be difficult to sys- 
tematically resum the large logarithms arising at higher 
orders from the phase-space restrictions. Therefore, to 
achieve the best theoretical precision, it is important to 
implement the central jet veto using an inclusive kine- 
matic variable. This allows us to derive a factorization 
theorem with analytically manageable ingredients, which 
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can then be used to sum large phase-space logarithms. 

We will consider a simple kinematic variable that ful- 
fills the above criteria, leaving the discussion of more so- 
phisticated generalizations to the next subsection. The 
key variables for the isolated Drell-Yan process arc shown 
in Fig. |31 The proton momenta Pg and P b are used to 
define lightlike vectors and n b , 

p£ = ^<> Pi l = ^<, as) 

where the protons are massless and n 2 = 0, n 2 = 0, 
and n a -n b = 2. Using the beam axis, we define two 
hemispheres a and b opposite to the incoming protons. 
We then divide up the total hadronic momentum as 

p1 c = BZ + BZ, (14) 

where B£ = p^ x and B b = p^ are the total final-state 
hadronic momenta in hemispheres a and b. Of these, we 
consider the components 

B+ = n a -B a = + tanht/ ) e~ 2y " , 

B+ = n b -B b = +tanhj/ b )e- 2 ^ , (15) 

where B® b are the energy components and y a>b are the to- 
tal rapidities of B^ b with respect to the forward direction 
n a _ b for each hemisphere. Here, lim y _ i . 00 (l +tanhy) = 2 
and 1+tanhy > 1.8 for y > 1, so B^ b scale exponentially 
with the rapidities y a ^ b . 

In terms of the measured particle momenta pk in hemi- 
sphere a, 

B a =J2 n "-P k = E^ 1 +tanh77 fe )e- 2 " fc . (16) 

Here, and rjk are the experimentally measured energy 
and pseudorapidity with respect to n a , and we neglect 
the masses of final-state hadrons. An analogous formula 
applies for B b . Hence, B+ and B b receive large con- 
tributions from energetic particles in the central region, 
while contributions from particles in the forward region 
are suppressed. Thus, requiring small B^ b <C Q is an 
effective way to restrict the energetic radiation in each 
hemisphere as a smooth function of rapidity, allowing 
forward jets and disallowing central jets. At the same 
time, soft radiation with energies <C Q is measured, but 
not tightly constrained. 

As an example, consider the cut 

Bl b <Qe' 2 y^. (17) 

This constraint vetoes any events with a combined energy 
deposit of more than Q/2 per hemisphere in the central 
rapidity region \y\ < y cut . In the smaller region \y\ < 
Vcut — 1, the energy allowed by Eq. (jTTJ) is reduced by a 
factor of e 2 ~ 7, essentially vetoing any jets there. In 
the larger region \y\ < y C ut + lj it is increased by the 
same factor, so beyond y cut + 1 the hadronic final state is 



hemisphere b 1 hemisphere a 




B£ = n b -b b + n b -k b i B+ = n a -b a +n a -k a 



FIG. 3: Definition of hemispheres and kinematic variables for 
isolated Drell-Yan. 

essentially unconstrained. Thus, a typical experimental 
value might be y cut = 2, which vetoes energetic jets in 
the central region \y\ < 1. The precise value of the cut 
on B^ b will of course depend on the requirements of the 
experimental analyses. 

Note that the variable is similar to the total trans- 
verse energy in hemisphere a, defined as 

E Ta = E — T— = E + tanh % )e"^ . (18) 
cosh?7fc ^-^ 

k£a k£a 

B+ has two advantages over Exa- First, the exponential 
sensitivity to rapidity is much stronger for which 
means it provides a stronger restriction on jets in the cen- 
tral region and at the same time is less sensitive to jets in 
the forward region. Second, since B£ is a specific four- 
momentum component and linear in four-momentum, 
(pi + p2) + = pi + P2 > h is much simpler to work with 
and to incorporate into the factorization theorem. It is 
clear that the isolated Drell-Yan factorization theorem 
discussed here can be extended to observablcs with other 
exponents, e~ ar)k , much like the angularity event shapes 
in e+e" [6J. 

One should ask, down to what values can B^ b be 
reliably measured experimentally? In principle, par- 
ticles at any rapidity contribute to B^ b , but the de- 
tectors only have coverage up to a maximum pseudo- 
rapidity ?7detj as indicated in Fig. [3] For the hadron 
calorimeters at the LHC ?7dct — 5 and at the Tevatron 
77<j e t — 4. In the hadronic center-of-mass frame, the 
unscattcred partons inside the proton have plus compo- 
nents of C(Aq CD /£' cm ), so any contributions from the 
unmeasured proton remnants are always negligible. The 
question then is, what is the maximal contribution to 
_Bj~, from initial-state radiation that is missed as it is 
outside the detector? In the extreme scenario where all 
proton energy is deposited right outside rjdet, we would 
have B+ b = UTeVe" 10 = 0.6 GcV at the LHC and 

B+ b = 2TcVe~ 8 = 0.7 GcV at the Tevatron. In more 
realistic scenarios, the contribution from such radiation 
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is suppressed by at least another factor of 10 or more. 
Therefore, the finite detector range is clearly not an is- 
sue for measuring values B^ b > 2 GeV, and the relevant 
limitation will be the experimental resolution in B^ b . 

The factorization theorem for isolated DrcU-Yan, 
which we prove in Sec. IIV1 reads 



1 



da 



cr dqMYdB+dB^ 



Y,H t] (q 2 ,fl) Jdk+dk+ (19) 



x q 2 Bi[u> a (B+ - fc+),x a ,^] J B i [w b ( J B 6 + - k+),x blf i] 



x Sn lcm i(ka , k^ , fi) 



1 + 



Aqcd UgjB+j 

Q 



Q 2 



The physical interpretation of Eq. (fT9|) is that we take 
partons i and j out of the initial-state jets Bi, Bj and 
hard-scatter them to final state particles with Hij , while 
including S'ihemi to describe the accompanying soft radi- 
ation. The hard function is identical to the one in 
the threshold factorization theorem in Eq. (TT2"j) , and the 
sum in Eq. (|19p is again only over ij = {uu, uu, dd, . . .}. 
The quark and antiquark beam functions B q and Bq de- 
scribe the effects of the incoming jets and have replaced 
the PDFs. The variables u> a> b = x a ^E cm . The hard par- 
tons are taken from initial-state jets rather than protons, 
so unlike in the threshold case the gluon PDF now con- 
tributes via the beam functions. We will see how this 
works in more detail in Sec. Ill CI Finally, S'ihemi is the 
initial-state hemisphere soft function. 

The kinematic variables in Eq. (fT9|) are displayed in 
Fig. [3] The soft function depends on the momenta 
k£ = n a ■ k a and k b — rib-kb of soft particles in hemi- 
spheres a and b, respectively. Much like PDFs, the beam 
functions Bi(t a ,x a , fi) and Bj(tb, x b , fi) depend on the 
momentum fractions x a and x of the active partons i 
and j participating in the hard collision. In addition, 
they depend on invariant-mass variables 



t a = uj a b+ > , t b = io b b£ > . 



(20) 



where uj a>b = x a j,E cm are the hard momentum compo- 
nents and b^ = n a -b a . The momentum 6^ is defined as 
the total momentum of the energetic particles radiated 
into hemisphere a, as shown in Fig. [3j and similarly for 
bf . (The kinematics are shown in more detail in Fig. 1131 ) 
Before the hard interaction, the momentum of the active 
quark can be written as 



(21) 



The first term is its hard momentum along the proton di- 
rection, and the last two terms are from the momentum 
it lost to radiation, where b 2 ± = —b 2 T contains the trans- 
verse components. The quark's spacclike invariant mass 



is —uj, 



ab + a -bl 



-t a 



b 2 T . The beam function Bi for 



x a E cm b+, which is 



hemisphere a depends on t a 
the negative of the quark's transverse virtuality. (When 
the distinction is unimportant we will usually refer to t 



simply as the quark's virtuality.) By momentum conser- 
vation fe+ = £>+ — fc+, leading to the convolution of the 
beam and soft functions as shown in Eq. (Tin)) . Physically, 
the reason we have to subtract the soft momentum from 
,B + is that the beam function only properly describes 
the collinear radiation, while the soft radiation must be 
described by the soft function. An analogous discussion 
applies to Bj and t b for hemisphere b. The convolutions 
in the factorization theorem thus encode the cross talk 
between the soft radiation and energetic collinear radia- 
tion from the beams. 

By measuring and constraining B+ we essentially mea- 
sure the virtuality of the hard parton in the initial state. 
As the proton cannot contain partons with virtualities 
larger than Ag CD , the initial state at that point must be 
described as an incoming jet containing the hard off-shell 
parton. This is the reason why beam functions describ- 
ing these initial-state jets must appear in Eq. (|19p . It 
also follows that since t 3> Aq CD we can calculate the 
beam functions perturbatively in terms of PDFs, which 
we discuss further in Sec. Ill CI 

It is convenient to consider a cumulant cross section, 
including all events with B^ b up to some specified value, 

as in Eq. (JTTJ) . Integrating Eq. (H]) over < B+ b < B+ ax 
we obtain 



1 der 
Vo dq 2 dY 

x Bi[u a [B\ 

X Sihemi(^|~! , fl) 



Aqcd 



1 + 



max 
Wa.bB 



(22) 



Q 



Q 2 



where the soft function S'ihemi is the same as in Eq. (jT9J) , 
and wc defined the integrated beam function 



Bi(t max , x,fi)= / dt Bi(t, x, n) 6»(t max - t) 



(23) 



The cut B+ b < B+ ax implies the limit b+ b < B+ ax - k+ b 
and t a ,b < w n ,i)(-Bmax — ^J' leading to the convolutions 
in Eq. (|22]l. 

The factorization theorem Eq. (|19|) and its integrated 
version Eq. ([22]) are valid in the limit t a ^/Q 2 — B^ b /Q = 
A 2 <C 1, and receive power corrections of 0(A 2 ). Thus, 
for i? max = Qe~ 2Vcllt with y cut = 1, we expect the power 
corrections not to exceed e~ 2 ~ 10%. This is not a fun- 
damental limitation, because the power corrections can 
be computed in SCET if necessary. If the soft function is 
treated purely perturbatively, there are additional power 
corrections of (Aqcd / B^ b ) , which account for soft sin- 
gularities as B^ b — > 0. The variables B\ are infrared 
safe with respect to collinear splittings [65{ . 

The hard function receives perturbative a s corrections 
at the hard scale fin — Q, the beam functions have 
a s corrections at the intermediate beam scale \j? b ~ 
t max ~ QB+ ax , and the soft function at the soft scale 
^5 ~ S+ ax . For example, for Q ~ 1 TcV and y cut = 2 
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we have /ib — 140 GeV and ^5 ~ 20GcV. Even with 
a very small Q ~ 100 GeV, perhaps for Higgs produc- 
tion, fiB — 14 GeV and [is — 2 GeV are still pertur- 
bative (although at this point nonperturbative contribu- 
tions ~ Aqcd / Ms to the soft function might no longer 
be small and may be incorporated with the methods in 
Refs. (ill, IE3])- I n fixed-order perturbation theory, the 
cross section contains large single and double logarithms, 
ln(B+ ax /Q) ~ -4 and \n 2 (B+ ax /Q) ~ 16, invalidating 
a hxed-order perturbative expansion. The factorization 
theorem allows us to systematically resum these loga- 
rithms to all orders in perturbation theory, which is dis- 
cussed in more detail in Sec. Ill Ft 




The factorization theorem Eq. ([19]) also applies to 
other non-hadronic final states such as Z' —> £ + £~, or 
Higgs production with H — > 77 or H — s- ZZ* —> 4£. In 
each case, q 2 and Y are the total non-hadronic invariant 
mass and rapidity, and central jets are vetoed with a cut 
on B^ b . The only dependence on the process is in the 
hard function, which must be replaced appropriately and 
can be taken directly from the corresponding threshold 
factorization theorem. One may also consider W produc- 
tion with W — ► Iv, with an appropriate replacement of 
q 2 and Y with the charged lepton's rapidity. For a light 
Higgs with Q ~ ma, the isolated Drcll-Yan factoriza- 
tion theorem applies to Higgs production through gluon 
fusion gg — > H and Higgs-strahlung qq — > VH, which 
are the dominant production channels at the LHC and 
Tevatron, respectively. 2 For a generic process pp — > XL, 
the sum over ij — {gg, uu, uu, dd, . . .} includes a gluon- 
gluon contribution, but still no cross terms between dif- 
ferent parton types, and there will be two independent 
soft functions S?^ emi and S^ cmi . [As shown in Sec. IIV1 
only the qq soft function contributes to isolated Drell- 
Yan, so the labels were omitted in Eq. Q19p.] Indeed, the 
gluon-gluon contribution involving the gluon beam and 
soft functions, B g and S?JL v gives the dominant contri- 
bution in the case of Higgs production. 

With the above physical picture, we can understand 
why the gluon beam function appeared in jp — > J/ipX 
in the analysis of Ref. [3(| in the limit where Ej/^ —> E~ r 
Taking px as the total momentum of final-state hadrons 
other than the J/ip, one has n-px ~ £"0111(1 — Ej/^/E^), 
where n is the proton direction. For Ej/^p close to E^, 
energetic radiation in the final state is restricted to a 
jet close to the n direction. Just as for our B~^ b , the 
measurement of Eji^ probes the radiation emitted by 
the colliding gluon in the initial state. Thus, the proton 
is broken apart prior to the hard collision, and the gluon 
beam function is required to describe the initial state. 



2 In vector-boson fusion and associated production gg — > ttH, the 
situation is more complicated and one has to explicitly consider 
the process pp — > XjjH with two forward (top) jets. 



FIG. 4: Generalized definition of hemispheres. The total ra- 
pidity of the leptons is Y, b^ b = n a ,b-b a ,b, and k^ h (Y) — 
n a ,b-k a ,b(Y). 



B. Generalized Observables 



The factorization theorem in Eq. (|19|) applies for t a <C 
q 2 and % <C q 2 ■ This includes the situation where in 
the hadronic center-of-mass frame there is a numerically 
significant asymmetry uj a = x a E cm > U)b = XbE CTa . 
This means that the boost between the hadronic and 
partonic ccntcr-of-mass frames, given by the leptonic 
Y = In ^/uja/uji, = In y/x a /xb, is significantly different 
from zero. We explore the implications of this here. 

If there is no hierarchy, uj a w w ,yu! a Wb — Q, cor- 
responding to Y ~ 0, we can define a simple variable to 
constrain both hemispheres simultaneously, 



B 

From Eq. (fl"9"|) . this gives 
1 da 



B+ + B2 



Q 



(24) 



dq 2 &YdB 



= ^H l3 {q 2 ,ll) fdt a dtb 
ij J 
X Bi(t a ,X a , n) Bj(t b ,X b , fl) 

xQSbIqB-^-^,^) , (25) 
V uj a uj b / 



where the soft function is defined as 
5s(fc + ,,u) = J dfc+dfc^ 5ihcmi(fca >^b ,n)d(k + -k+-k 



(26) 

The advantage of using B is that the soft function now 
only depends on the single variable k + = k+ + k b , much 
like the soft function for thrust in e + e~ collisions. 

If we have a hierarchy ujb < Q < 0J a , the final state has 
a substantial boost in the n a direction, as shown in Fig. [4] 
In this case, the energetic radiation will generically be 
much closer to the beam axis in hemisphere a than in 
hemisphere b. To take this into account, it is natural to 
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impose different cuts on £?+ and Bt . Using the boost- 
invariant combinations ui a B+/q 2 and UbB b /q 2 to define 
the cut, we obtain 



R+ 

= —2- < e - 2 v<= 

UJ b ~ 



UbBt B b 



< e 



-2y cut 



(27) 

so B+ has a tighter constraint than B b , as desired. If we 
simply replace B by 2?+ /ujb + B b / ui a , the soft function 
analogous to Sb in Eq. (|26p will depend on the combina- 
tion (uj a k+ + uj b k^)/Q 2 . 

However, we should also adjust the hemispheres them- 
selves to take into account the significant boost of the 
partonic center-of-mass frame. We therefore define a 
generalized hemisphere a as y > Y and hemisphere b 
as y < Y, as shown in Fig. [U The corresponding total 
hemisphere momenta are denoted as B* b (Y) and the soft 
hemisphere momenta as k^ b (Y). The original definitions 
in Fig.[Ucorrespond to B+ b (Q) = B+ b and k+ b (0) = k+ b . 

The generalization of B is given by the boost-invariant 
combination 



tb = 



cj a B+(Y)+u b B+(Y) 



r 



(28) 



With the generalized definition of the hemispheres, 
B^ b (Y) and ui a .b transform under a boost by y in the 
n„ direction as 



Bt{Y) 



U>b 



BfiY + y) 
ui a = e y u} a , 



e- y B+(Y) 
e y B+(Y), 



i. — f> y 



UJb 



(29) 



Thus, boosting by y = —Y from the hadronic to the 
partonic center-of-mass frame gives 

B+'W + B+'jo) r „ m 

tb = 5 = -p. • (oU) 



In the partonic center-of-mass frame we have ui' a = ui' b = 
Q, so there is no hierarchy. Correspondingly, the gener- 
alized hemispheres in this frame are again perpendicular 
to the beam axis, so Eq. (f5U|) has the same form as B. 

Note that for e + e~ — > jets, one can use the thrust 
axis to define two hemispheres with n a ,b analogous to our 
case. In the 2-jet limit, thrust is then given by 1 — T = 
(Qn a -px a + Q n b -px b )/2Q 2 - Hence, we can think of tb 
as the analog of thrust for incoming jets. For this reason 
we will call tb the "beam thrust" . 

In analogy to Eqs. (jTTJ) and ([27l) . we define the cutoff 
on t b by 



allow energetic particles up to rapidities y < y™* (with 
y measured in the partonic center-of-mass frame). 

The beam functions are boost-invariant along the 
beam axis, so the different hemisphere definitions 
do not affect them. The soft function is boost- 
invariant up to the hemisphere definition, which de- 
fines its arguments k^ b . Hence, boosting by —Y 



we have S , ihcmi [e r fc+, e Y k+;Y] = 5 ihcmi [fc+, k+ ; 0] = 
5ihcmi(fc^ , fcjj"), where the third argument denotes the 
definition of the hemispheres. This implies that the soft 
function for tb is the same as in Eq. (|26p . The factoriza- 
tion theorem for tb following from Eq. (|19p is 



1 



da 



ctq dq 2 dYdTB 



ij J 

x Bi(t a , x a , /i) Bj(tb, Xb, /i) 
x Q S b (Qt b - — ~t tb , n 



Integrating over < tb < exp(— 2y B nt ) we obtain 



dq 2 dY [VB 



e X p(-2 lJ ™ t ) 



dr 



da 



D 



dq 2 dYdr B 



(32) 



(33) 



We will use Eqs. (|32j) and (|33|) to show plots of our results 
in Sec. E3 



C. Relating Beam Functions and PDFs 

The beam functions can be related to the PDFs by per- 
forming an operator product expansion, because t a ,b 3> 
Aq Cd . 3 This yields the factorization formula 



Bi(t,x,fx) = 



i + o 



A 



QCD 



(34) 



where we sum over partons j = {g, u, u,d, . . .}, Iy are 
perturbatively calculable Wilson coefficients, and fj is 
the standard PDF for parton j. The £>(A| CD /t) power 
corrections in Eq. (|34p involve proton structure functions 
at subleading twist. Further mathematical details on 
Eq. ([34)) are discussed in Sec. Mil whereas here we fo- 
cus on the physical ramifications. 

The interpretation of Eq. ([34"]) is illustrated in Fig. [5] 
At a hadronic scale fj,\ ~ 1 GeV, the initial conditions for 
the PDFs fj can be specified, and one has the standard 
DGLAP evolution up to the scale /x^, 



(35) 



tb < e 



(31) 



For tb — >• or equivalently y B ut — > oo the jets along the 
beam axes become pencil-like, while for generic y c B nt we 



A detailed discussion of the appropriate operator product expan- 
sion is given in Ref. |50|| . 
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fl\ changing x fig changing t flfj 

FIG. 5: Evolution of the initial state. Starting from the low scale [ia, the incoming proton is described by the ir-dependent 
evolution of the PDFs, which redistributes the total momentum of the proton between its constituents. At the scale /is, the 
proton is probed by measuring the radiation in the final state and breaks apart. This is the scale where the PDFs are evaluated 
and the a;-dependent evolution stops. Above /is, the proton has ceased to exist, and the initial state behaves like an incoming 
jet, whose evolution is governed by the virtuality t of the off-shell spacelike parton that eventually enters the hard interaction 
at the scale /j,h- 



The anomalous dimensions Pjy are the standard QCD 
splitting functions for quarks, antiquarks, and gluons (in- 
cluding the color factors and coupling constant). Equa- 
tion ([34]) applies at the scale fi = hb, since this is the 
scale at which a measurement on the proton is performed 
by observing the soft and collinear radiation contribut- 
ing to B^ b . At this scale, a parton j with momentum 
fraction £ is taken out of the incoming proton according 
to the probability distribution /j(£,/i). As the parton 
continues to propagate and evolve with pi > fis, it is 
modified by virtual radiation and by the emission of real 
radiation, which forms a jet. The evolution in this region 
no longer depends on £, but instead on the virtuality t. 
This evolution occurs with fixed x and fixed parton type 
i, via the beam function RGE 

d f 

fi—B l {t,x,n)= J dt'y B {t-t',n)B t (t',x 7 fi). (36) 

This result for initial-state jet evolution has the same 
structure as the evolution for final-state jets. In fact, the 
anomalous dimension 7^ is identical to that for the quark 
jet function to all orders in perturbation theory [Ho[. We 
discuss this correspondence further in Sec. Mil 

The effect of initial-state real and virtual radiation is 
described by the perturbatively calculable Wilson coef- 
ficients Iij(t, x/£, /i) at the scale \x = \ib- They encode 
several physical effects. The virtual loop corrections con- 
tribute to the Xu and modify the effective strength of the 
various partons. If the radiation is real, it has physical 
timelikc momentum. Hence, it pushes the active parton 
in the jet off shell with spacelike virtuality — t < and 
reduces its light-cone momentum fraction from £ to x. 

In addition, the real radiation can change the iden- 
tity of the colliding parton, giving rise to the sum over 
j in Eq. (|34[) . For example, an incoming quark can ra- 
diate an energetic gluon which enters the hard interac- 
tion, while the quark itself goes into the final state. This 



gives a contribution of the quark PDF to the gluon beam 
function through X gq . Similarly, an incoming gluon can 
pair-produce, with the quark participating in the hard 
interaction and the antiquark going into the final state. 
This gives a contribution of the gluon PDF to the quark 
beam function through X qg . There are also of course real 
radiation contributions to the diagonal terms, T qq and 
1 gg , where the parton in the PDF and the parton partic- 
ipating in the hard interaction have the same identity. 

At lowest order in perturbation theory, the parton 
taken out of the proton directly enters the hard inter- 
action without emitting radiation, 

4 ree (*>!^) =M(*)*(l"|)- (37) 

Thus at tree level, the beam function reduces to the PDF 

B?~(t,x,ri = 5(t)f i (x,ii). (38) 

Beyond tree level, Tij(t,x/£, /i) can be determined per- 
turbatively as discussed in more detail in Sec. IIIH where 
we give precise field-theoretic definitions of the beam 
functions and quote the one- loop results for X qq and T qg . 

Interestingly, in the threshold factorization theorem 
Eq. (fT2|) . cross terms between quark and gluon PDFs 
are power suppressed, so the gluon PDF does not con- 
tribute at leading order. In the inclusive case Eq. (J9)), 
such cross terms are leading order in the power counting. 
For isolated Drell-Yan, there are no cross terms between 
quark and gluon beam functions, but there are leading- 
order cross terms between different PDFs, which appear 
via the contributions of different PDFs to a given beam 
function in Eq. (f3"4")) . Thus, the isolated case is again 
in-between the inclusive and threshold cases. 
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D. Comparison with Initial-State Parton Shower 

The physical situation associated with the beam evo- 
lution has an interesting correspondence with that of 
initial-state parton showers. As pictured in the region 
between /xb and \ih in Fig. [5l the parton in the beam 
function evolves forward in time while emitting a shower 
of radiation into the final state governed by the anoma- 
lous dimension j l B (t — t' , /i) in Eq. (|36[) . This equation 
has no parton mixing. Each emission by the radiating 
parton increases the magnitude of its spacelike virtuality 
—t < 0, pushing it further off-shell in a spacelike direc- 
tion. At the time the parton is annihilated in the hard 
collision, it has evolved to some t with \t\ q 2 , so the 
large momentum transfer q 2 guarantees that no partons 
in the final state are spacelike. This description agrees 
quite well with the physical picture associated with the 
evolution of the primary parton in an initial-state parton 
shower, as summarized in Ref. |36| . 

Differences in the description arise when one considers 
the initial-state parton shower in more detail (for sim- 
plicity we focus on the so-called longitudinal evolution). 
The shower is based on the evolution equation for the 
PDFs in Eq. (|3"5j) . An evolution forward in time is not 
practical because of the lack of prior knowledge of the 
scale of the hard interaction, so the shower uses backward 
evolution starting at a given partonic hard scale Q [47j |. 
Knowing the identity of the final parton i, the shower 
evolves based on the probability dPi/dt that parton i is 
unresolved into parton j via the splitting j — > if at an 
earlier (lower) scale t. The evolution equation is (36j 



dT^i (x, tmaxi ^) 

dt 



E 

jj' 

x -Viix,^ 



P,-. 



,(z,t) 



ft(x,t) 



(39) 



where Vi(x,t ma , x ,t) is the shower Sudakov exponential, 
which is interpreted as the probability for no emissions to 
occur between the initial value t max and t. The evolution 
variable t, which determines the scale of the splitting, is 
usually chosen as the virtuality or transverse momentum 
of the parton. 

The mixing of partons in the PDF evolution influences 
the shower. In particular, the evolution kernel depends 
on the PDF fj (x/z,t), which determines the number den- 
sity of partons of type j at the scale i, and inversely 
on the PDF fi(x,t). Thus, unlike in the beam evolu- 
tion in Eq. (|36|) . the shower evolution in Eq. ([39| still 
knows the identity of the initial-state hadron. Double 
logarithms in the initial-state parton shower are gener- 
ated in q — > qg and g — > gg splittings because of the 
soft-gluon singularity ~ 1/(1 — z) in the splitting func- 
tions. This singularity is regulated [36| by the upper 
cutoff z max = xj(x + x e ), where x e provides a lower cut- 
off on the gluon energy in the rest frame of the hard 
scattering, E g > x e -fE cnl /2 ~ 2 GcV (where 7 is the 
boost factor of the hard scattering). Hence, one loga- 



rithm, lnx £ , is generated by the z integration, and one 
logarithm, lni, by the collinear 1/t singularity. In con- 
trast, the beam function contains double logarithms In 2 t 
similar to a final-state parton shower, where the z inte- 
gration yields a kernel ~ (In t) /t that produces a double 
logarithm In 2 t via the t evolution. 

The above comparison is very rough. For example, 
the influence of soft radiation on both the shower and 
on the isolated factorization theorem was not compared 
and is likely to be important. Furthermore, the goal of 
the shower is to provide a universal method for populat- 
ing fully exclusive final states, while the beam function 
applies for a more inclusive situation with a particular 
measurement. Note that just the presence of mixing in 
the initial-state parton shower and absence of mixing in 
the beam-function evolution does not imply an inconsis- 
tency. For example, it is well known that the final-state 
parton shower reproduces the correct double logarithms 
for e + e~ event shapes [HI, even though there is no par- 
ton mixing in the evolution of the corresponding hard, 
jet, and soft functions. In the future it would be inter- 
esting to test in detail the correspondence between the 
double logarithms generated by the initial-state parton 
shower and those predicted by our factorization theorem 
for the isolated Drell-Yan process. 



E. Relation to Fixed-Order Calculation 

The factorization theorem for the cross section in 
Eq. (|19|) and the factorization for the beam function in 
Eq. (|34jl together allow us to describe in more detail how 
various Feynman diagrams that would appear in a fixed- 
order calculation contribute to the cross section in our 
kinematic region. Various examples are shown in Fig. [6] 
In Fig. 6(a) we have the tree-level qq annihilation pro- 
ducing a 7 or Z, which involves the tree-level 0{oP s ) hard 
function, beam functions, and soft function, denoted by 
a superscript (0) in the figure. In Fig. |6(b)[ initial-state 
gluons couple to a quark loop (e.g. a top quark), which 
subsequently annihilates into a 7, Z, or Higgs. The 
quarks in this loop are far off shell, so they can be inte- 

to 



grated out and appear as one-loop corrections, H, 
the hard coefficient in the factorization theorem. Other 
possibilities for this graph are power suppressed. 

The situation for the vertex correction in Fig. 6(c) is 
more involved. If the gluon in the loop is hard, all par- 
ticles in the loop are far off shell and can be integrated 
out, giving the one-loop hard function shown as the 
first term on the right-hand side. In the second term, the 
gluon is collinear to the incoming quark beam and gives 
a virtual one-loop contribution to the quark beam func- 
tion, Bq . The third term is the analog of the second, 
but now with the gluon collinear to the incoming anti- 
quark. Finally in the fourth term, the gluon is soft, com- 
municating between the incoming collinear beams. Here, 
the eikonal approximation holds for describing the quark 
propagators. The generalization of this to all orders in 
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FIG. 6: Factorization for isolated Drell-Yan in pictures. The left-hand side of each equality are graphs in QCD, while the 
right-hand side shows the sum of the corresponding SCET diagrams. Dashed lines are collinear quarks, and springs with a line 
through them are collinear gluons. The double lines denote soft Wilson lines, and the gluons attached to them are soft. 



a s leads to the fact that the soft function is a matrix 
element of Wilson lines. Although a single loop graph 
contributes in several different places in the factorization 
theorem, all of these contributions have a precise separa- 
tion in SCET. We will use this separation in Sec. IIVI to 
prove the isolated Drell-Yan factorization theorem. 

An interesting contribution occurs in Fig. |6(d)| where a 
gluon is radiated into the final state. Because of the kine- 
matic restrictions in isolated Drell-Yan, this gluon can 
only be collinear to the incoming quark, collinear to the 
incoming antiquark, or soft, and these three possibilities 



are represented by the diagrams on the right-hand side of 
the equality. In the first case, we have a real-emission cor- 
rection to the quark beam function, Bq X \ In the second 
case, the intermediate quark is far off shell and can be 
integrated out, and the gluon collinear to the antiquark 
arises from a collinear Wilson line contribution in Bq~\ 
The third case gives a real-emission correction to the soft 
function, sj^ . The full-theory graph in Fig. |6(dj] has a 
t-channel singularity. An important fact about the iso- 
lated Drell-Yan factorization theorem is that it fully cap- 



13 



tures the dominant parts of this singularity, and allows a 
simple framework for a rcsummation of higher order a s 
corrections enhanced by large double logarithms due to 
this singularity. For threshold Drell-Yan, the kinematic 
restrictions are stronger and only allow the third graph 
with soft initial-state radiation. In inclusive Drell-Yan, 
the gluon is treated as hard, and the graph in Fig. |6(d)| 
only corrects , without providing a framework for 
summing the large double logarithms that appear when 
we make a global measurement of the radiation in each 
hemisphere defined by the beams. 

The situation is a bit simpler for Figs. 6(e) and |6(f)| 
In Fig. 6(e) the incoming collinear gluon from the PDF 



pair-produces a quark and antiquark both collinear to 
this beam direction, and the quark enters the hard inter- 
action. Therefore, this is a one-loop correction to the 
quark beam function, B q x \ proportional to the gluon 
PDF f g . The beam functions again allow us to resum 
the possibly large logarithms due to this f-channcl singu- 
larity. Other possibilities for the final-state antiquark in 
Fig. 6(e) lead to power-suppressed contributions. Sim- 



ilarly, the s-channel graph in Fig. 6(f) which has the 
same initial and final states as Fig. 6(e) has no leading- 



power contribution and only contributes to Eq. (|19p in 
the power-suppressed terms. The same is also true for 
Drell-Yan in the threshold region. Only inclusive Drell- 
Yan receives a leading-order hard contribution from the 
s-channel graph, which is then treated as of the same size 
as the i-channcl graphs. 



F. Renormalization Group Evolution 

In this subsection, we discuss and compare the struc- 
ture of large logarithms in the cross sections for inclu- 
sive, threshold, and isolated Drell-Yan. These large loga- 
rithms may be summed using the renormalization group 
evolution of the individual functions appearing in the fac- 
torization theorems. In fact, the structure of large log- 
arithms in the differential B* b cross section allows us 
to infer the necessity of the beam functions in the iso- 
lated factorization theorem. This procedure provides a 
method of determining whether beam functions enter for 
other observables or processes than those studied here. 
The consistency of the RGE was used to provide a sim- 
ilar consistency check in Ref. [69| when deriving a new 
factorization theorem for the invariant-mass distribution 
of jets initiated by a massive quark in e + e _ collisions. 
In that case, the RGE consistency provided important 
constraints on the structure of the factorization theorem 
at scales below the heavy-quark mass. 

In inclusive Drell-Yan, the hard functions -ff]" cl are sen- 
sitive to the scale ~ Q of the hard interaction, and 
the proton mass defines a low scale ha — 1 GeV > Aqcd 
(which is still large enough so perturbation theory can 
be applied for the PDF evolution). The measurement 
of q 2 and Y in this case does not introduce additional 
scales, and thus does not influence the structure of the 



logarithms. Thus, we have the hierarchy ha <C Hh, and 
the large logarithms are L = \yi(ha/ Hh)- Here, only 
single-logarithmic series, (a s L) k , arc generated at higher 
orders in perturbation theory. The logarithms are factor- 
ized as ln(/// hh) + ^-(ha/ h) in the factorization theorem 
in Eq. ^ and may then be resummed. The general form 
of the running is pictured in Fig. [7J a ) ■ The logarithms 
ln(/ZA/V) are summed by evolving the PDFs /i(£ OJ h) and 
fj{£b,n) from ha up to the common scale fi- The inclu- 
sive hard function, H lncl (x a /^ a , Xb/^b, Q 2 , A 4 ), is evolved 
from hh down to fx, summing the logarithms ln(/// Hh)- 
The choice of h is arbitrary. Taking /i ~ hh corresponds 
to only running the PDFs up, while for (j, — V>A only H lncl 
runs down. The equivalence of these two choices implies 
that H lncl must be convoluted with the two PDFs and 
exhibit a factorized structure for logarithms in the a and 
b variables. 

For threshold Drell-Yan, the kinematic restrictions 
only allow soft radiation in the final state. This induces 
additional large logarithms ln(l — r). These can be writ- 
ten in terms of a ratio of scales ln(/xs//x#), where the 
soft scale hs — Q(l — t) is another important scale in 
the analysis. The logarithms L = \tl{hs / Hh) appear as 
double-logarithmic series (a s L 2 ) k in the cross section. 
In the threshold factorization theorem in Eq. (|12p . these 
double logarithms can be summed by evolving the PDFs 
and the threshold soft and hard functions, Sthr and H, to 
a common scale Hi as shown in Fig. Efb). Since ^ a ,b — > 1, 
the logarithms ln(l — £ a ) and ln(l — £&) are also large. 
The RGE for the PDFs must be expanded, and the result 
sums a double-logarithmic series of ln 2 (l-£) terms. The 
threshold soft function sums double logarithms ln 2 (/i//is) 
between hs and Hi while the threshold hard function 
sums double logarithms \t?{hI Hh) between hh and \x. 
The evolution equations are 



H—H(q 2 ,H) =-fH{q 2 ,H)H(q 2 ,H), 

"H£ /i( *'" ) = /y^ xpanded (|"")/^'"")' ( 4 °) 

The consistency of the RGE at the scale \i shown in 
Fig. [T^b) implies that the double logarithms in /j, fj, 
and Sthr combine in such a way that the RGE of the 
convolution fifj <S> Sthr is identical to that of H, and 
hence only depends on q 2 . 

For isolated Drell-Yan, the kinematic restrictions al- 
low both soft and collinear initial-state radiation, and 
induce an invariant-mass scale for each beam function, 
H 2 B — x a E cm B+ and n 2 B ~ XbE cln B^ , and a soft scale 
Hs — B+ b - For simplicity, we use a common scale hb 
for both beam functions in our discussion here. (Since 
the evolution of the two beam functions is indepen- 
dent, one can just as easily implement two indepen- 
dent beam scales.) As we saw in Sec. Ill Al at partonic 
center-of-mass energies of a hundred GeV to a few TeV 
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FIG. 7: RGE running for different Drell-Yan scenarios. Case (a) corresponds to the inclusive case. Case (b) corresponds to the 
threshold case, where the kinematics forces all hadrons in the final state to be soft. Case (c) corresponds to the isolated case. 
Here, the PDFs freeze out at the intermediate beam scale ^b, above which they are replaced by beam functions. 



there is a large hierarchy between the different scales, 
MA <C /is <C I^b ^ H>h, and correspondingly large dou- 
ble and single logarithms of the ratios of these scales. 
The RGE running for this case is shown in Fig. [7]Jc). 
Here, the PDFs are not restricted to their endpoints, so 
their evolution is given by Eq. (|35|) . which involves the 
uncxpanded and nondiagonal Pi j {£,/£!) and sums single 
logarithms, (a s L) k . For each fj this evolution joins at 
fi = fiB with the Wilson coefficients 2^ in the beam 



double logarithms in H . Thus, the structure of double 
logarithms necessitates the presence of beam functions in 
the isolated factorization theorem. 

By the same argument we can conclude that for all 
processes involving a threshold-type hard function H 
with double logarithms, and with x a ,b away from one, 
the description of the initial-state radiation will require 
beam functions B;. This includes all situations where H 
is the square of Wilson coefficients of SCET operators, 



function factorization B l = lij <g> fj of Eq. ((M]). The I,j H = J2k\Ck\ 2 (for example when the energetic partons in 

the hard collision all have distinct collinear directions). 
In particular, the theoretical description of any threshold 
process with x — > 1 can be extended to a factorization 
theorem for the respective isolated case with x away from 
one. This is achieved by adding variables B^ b , replacing 
the PDFs by beam functions, and replacing the thresh- 
old soft function by an appropriate soft function for the 
isolated case. 

Thus, beam functions are quite prevalent for cross sec- 
tions that one may wish to study at the LHC. In situa- 
tions where the hadronic final state is constrained with 
variables that are more complicated than B^ b , one gener- 
ically expects to find different beam functions and differ- 
ent soft functions encoding these constraints. This ex- 
tension is analogous to how the choice of jet algorithm 
modifies the definition of the jet and soft functions for 
central jets produced by the hard collision [17|. Even 
with this generalization, the beam and soft functions will 
both sum double-logarithmic series, and we expect that 
the factorization relating the beam function to the PDFs 
will carry through, just with different coefficients 1^. 



cancel the ^-dependent evolution of fj, and turn it into 
the i-dependent evolution of B i7 which sums a double- 
logarithmic series. The objects meeting at the common 
scale /i in Fig. (TJ^c) are the hard function, which is iden- 
tical to the threshold case in Eq. (f4T)|) , and the beam and 
soft functions, 

d /* 

H—B^t^x^) = J At' ^(t - t 1 \n) Bi(t',x,fj,) , 

M j^'5'ihcmi(fc c ^ I j M) = J "dfc Q dk b S[i lcm i(k a , k bl /i) 

x 1S iitetai (ki - K,^ - k' b , fi) . 

(41) 

The consistency of the RGE at /i now implies that the 
double-logarithmic running in the different variables for 
Bi, Bj, and Sihemi cancels such that the convolution 
BiBj (g) Sihemi has an RGE identical to H, which only 
depends on q 2 . (A detailed discussion of this consistency 
can be found in Ref. [7(| for the analogous case of two 
jet functions and the final-state hemisphere soft function, 
J JgifShemi, and in Ref. [5(| for the case discussed here.) It 
is important that this cancellation would not be possible 
if we tried to replace Bj by fi in the isolated factoriza- 
tion theorem. Given the type of double logarithms in 
the cross section, the single logarithms summed by the 
PDFs at generic x cannot combine with the double log- 
arithms in Sihemi to give a result in agreement with the 



G. Extension to Final-State Jets 

In Fig. [21 where we show the types of hadronic final 
states for inclusive, endpoint, and isolated Drell-Yan, 
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FIG. 8: RGE running for dijet production for (a) the threshold situation and (b) the isolated situation. 



we also included analogs where the lepton pair is re- 
placed by two jets. Figure |2(d)| shows the threshold 
dijet production process studied in Ref. [§1, which is 
the generalization of the Drcll-Yan threshold process in 



Fig. 2(b) Figure 2(e) shows the isolated dijet production 



process, which is the generalization of isolated Drell-Yan 
in Fig. 2(c) The goal of this subsection is to give a rough 



■^-variable of Drell-Yan are 



M 2 JJ = (p 1 + p 2 y 



2Pi • P 2 = M]j - M[ - M\ 



(42) 



where M jj is the total invariant mass of the two jets, and 
m 2 jj is their total invariant mass minus their individual 
invariant masses. The corresponding analogs of the Drcll- 



Here <f> and rj denote the azimuthal angle and pseudorapidity. Al- 
though our notation corresponds to a cone algorithm, at the level 
of our discussion we may equally well substitute a fcy algorithm. 
We take P± 2 t° be the momenta of the two hardest jets found by 
the jet algorithm. We will see below that the final-state restric- 
tions considered here eliminate the possibility of having other 
hard central jets. 



Yan t- variable are then 



tj 



M 2 
E 2 



At,/ = 



•jj 



E 2 



(43) 



idea of how the isolated factorization theorem will look 
when it is extended to include final-state jets. Recall that 
our proof of factorization is only for the Drell-Yan case. 
We stress that the factorization formula for the isolated 
dijet case discussed here expresses our expectations and 
has not been rigorously derived. 

Final-state jets are identified by a jet algorithm as 
more-or-less isolated groups of energetic particles within 
a cone 4 of some radius R = [{Acf)) 2 + (Arj) 2 } 1 / 2 . For a j^ e £ || 6,7 
dijet event, the jet algorithm allows us to define the total 
jet momenta Pf and P 2 - 5 Given these, we let y\ and 
i/2 be the rapidities of the two jets relative to the beam 
axis and define Ay = y\ — y 2 . The invariant masses of 
the jets are denoted as M 2 2 = P\ 2 - Two analogs for the 



Identifying the two jets already restricts the hadronic 
final state, which means there is no analog of the inclusive 
Drell-Yan factorization theorem. In the threshold case, 
we take the limit tj — > 1, which ensures that the final 
state consists of two back-to-back jets plus soft radiation 



and no additional energetic jets, as shown in Fig. 2(d) 
The limit Atj — > 1 is even more restrictive, since it also 
forces the two jets to have very small invariant masses, 
essentially behaving like massless particles for the factor- 
ization theorem. 

The threshold factorization for rrijj was considered in 



da 



^H IL (m 2 jj, Ay) [ d£ a d& dpi d P2 

cm J 



dm 2 JJ d(Ay) El 
x /(£a)/(&>) Jthr(mjj,pi,R) Jthi(mjj,p 2 ,R) 
x Sftte [mjj(l- Atj- (1-&) - (1-&) - Pl -p 2 ), Ay] 



1 + o f±0£E ,R,1-Arj 



(44) 



where one sums over the color basis IL, and for sim- 
plicity the dependence on flavor labels and p of the var- 



To the best of our knowledge a proof of the decoupling of Glauber 
gluons does not exist for threshold dijet production in hadronic 
collisions, so there is no complete proof of Eq. 1441 . 
When replacing m?jj by Mjj in Eq. 14411 , the main difference is 
that now «; = 2k°/Mjj. In this case, the threshold limit tj — > 1 
alone does not constrain the jet invariant masses, M?, to be 
small. Since M? ~ R 2 Mjj, they are constrained to be small by 
jet algorithms with R 2 <C 1. This induces complications in deriv- 
ing an all-orders factorization theorem, but still suffices to imply 
that the factorization formula in Eq. [144 t i with the replacement 
will sum the next-to-leading logarithms [f| . 
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ious functions have been suppressed. 8 The first argu- 
ment of the soft function 5|j{ hr is the energy it radiates 
outside the jet cones. The jet function J tri r depends on 
Pi = 2k®/mjj + Mf/rrijj, where k° is the energy of par- 
ticles it radiates outside its cone. The threshold limit 
At,/ — > 1 forces £ a> b — > 1 and pi i2 — > 0, so we have 
PDFs, jet functions, and soft functions that all corre- 
spond to the threshold limit and contain double loga- 
rithms. 9 When these functions are convoluted they con- 
sistently reproduce the double logarithms encoded in the 
renormalization group evolution of H IL , which is shown 
in Fig. [5]Ja). Compared to threshold Drell-Yan as de- 
scribed in Sec. Ill F[ there are extra convolutions for the 
jets and a more sophisticated soft function that is a ma- 
trix in color space. 

To extend the isolated Drell-Yan factorization theorem 
to the dijet case, we need to define analogs of the B^ b 
variables in Sec. Ill Al that can constrain the final state in 
an appropriate manner. We first define lightlike vectors 
along each jet direction, = (l,ni) and n 2 = (1,??2) 
with n 12 = Pi,2/\Pi,2\, and corresponding lightlike vec- 
tors in the opposite directions, n^ 2 = (1, — ^1.2)- Next, 
the final-state particles arc divided into four categories 
{a, 6, 1,2} as follows. We define TZi and 72-2 as all par- 
ticles that have been grouped into jets 1 and 2 by the 
jet algorithm. The remaining particles not grouped into 
either of 1Z 1,2 are divided into the two hemispheres a 
and b as before, which defines !Z a ,b- We then define plus 
momenta 



Pi = ^2 ni -p k , P 2 = Y2 712 ' Pk ' 

B+ = £ n a - Pk , B+ = n b-Vk ■ (45) 



k£ll a 



ken b 



This definition of B^ b is identical to Eq. ([2]). Since the 
jet algorithm is used for the grouping, these variables are 
infrared safe. The union of the four categories covers all 
of phase space, so the measurement of all the momenta 
in Eq. (|45jl defines a global observable sensitive to all ra- 
diation in the event. Just as in our Drell-Yan discussion, 
the definition of the variables B^ b ensures that radiation 
outside the reach of the detector can safely be ignored. 

For the isolated dijet limit we demand that 
Pf 2 /Mjj < 1 and B+ b /Mjj < 1. In addition, we con- 
strain y\ and y 2 such that the jets lie in the central re- 
gion sufficiently separated from the beam directions. The 
condition on P x + 2 ensures that the jet regions IZ1.2 only 
contain energetic radiation along the direction of their 



8 We also made a redefinition so that the PDFs / depend on light- 
cone momentum fractions rather than fixed energy as in Ref . Q , 
absorbing the difference into the hard functions Hjl . 

9 We follow the SCET definition of soft functions as matrix ele- 
ments of eikonal Wilson lines without subtractions, so S^jthr nas 
double logarithms. In SCET the jet functions have subtractions. 



jet plus soft radiation. The condition on B~f b has a sim- 
ilar effect as for isolated Drell-Yan. It ensures that there 
is only soft and no energetic radiation in the central re- 
gion apart from the two jets. Thus, we have exactly two 
isolated central jets. 

Since each category {a, b, 1, 2} predominantly contains 
all the corresponding collinear particles, this division of 
phase space mainly affects how the soft radiation is asso- 
ciated to each jet. In analogy to isolated Drell-Yan, we 
divide the total soft momentum as k = k a + fc& + k\ + fc 2 , 
where each fej is the total momentum of soft particles in 
TZi , and we define kf = n^ki . The corresponding isolated 
dijet soft function, S^(k^, k b , k^, k 2 , yi, 2/2)1 depends 
on all four directions ni, n 2 , n a , nb and hence on the 
rapidities y\ and y 2 . It now contains both incoming and 
outgoing soft Wilson lines and is a matrix in color space, 
where we can use the same color basis {LI} as in the 
threshold case. The soft function itself of course differs 
from the threshold case. 

The total jet momenta can now be written as 



1 2 



2 



(46) 



and similarly for P 2 . The first three terms on the right- 
hand side are the contributions from the energetic radi- 
ation in the jet, with wi, 2 ~ Mjj and q^,q± <C Mjj. 
Expanding in the small components of the beam and jet 
momenta, the hard momentum components in the beams 
and jets have to satisfy 



uj b - 



LOl- 



UJ2- 



(47) 



Thus, in this limit the two jets are massless and back- 
to-back in the transverse plane, but need not be back-to- 
back in three dimensions. In terms of ?/i i2 and Mjj = 
uj a LOb — LO\LOin\ -712/2, we then have 

uj a = Mjje^ + y^/ 2 = x a E cm , 
io b = Mjje-^ + y^/ 2 = x h E cm , 
cosh yi_ 2 



Wl,2 



M 



jj 



cosh(Ay/2) 



(48) 



The collinear radiation in the jets is described by jet 
functions that depend on the invariant-mass variables 
oj\ql and u>2q 2 ■ By momentum conservation we have 
q~l 2 = — ^1^2: so the jet functions are convoluted 
with the soft function through kf 2 . The subtraction of 
lS necessary to remove the plus momentum of soft 
particles in the jet, since the momentum distribution of 
these particles is properly described by the soft function 
not by the jet function. Just as for isolated Drell-Yan, the 
collinear initial-state radiation in the beams is described 
by beam functions, which depend on the invariant-mass 
variable t a = uj a {B+ — fc+) and momentum fraction x a , 
and similarly for hemisphere b. Since the jets are well- 
separated from the beams, removing the particles in 7Z± t 2 
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from the hemispheres H a ,b mainly affects the soft radia- 
tion and not the energetic partons collincar to the beams. 
Therefore up to power corrections, we expect the same 
inclusive beam functions as before. 

From the above discussion it is natural to suppose that 
the factorization theorem for isolated dijet production 
will be 



da 



dMjj Ay i Ay 2 AB+ AB+ AP+ AP+ 



(49) 



H IL (Mjj, y u y 2 )J dfc+ dk+ Ak+ Ak+ 

x J<MP+ - k+)] J<[wa(P a + - k+)} 

x B[uj a (B+ - k+), x a ] B[u b (B+ - k+), x b ] 

x S^ Ts (k+,k+,k+,k+, yi ,y 2 ) 

Aqcd 



1 + 



Ua,bB+ 



M 



jj 



M 2 
iv1 jj 



M 2 
m jJ 



where we again suppressed flavor labels and p, depen- 
dence. The hard function H IL (Mjj,y{) is precisely the 
threshold hard function, and we sum over the same color 
basis {IL}. The subscript < on the soft and jet functions 
denotes the fact that their plus momenta depend on the 
regions TZi, which in turn depend on y^. 

The consistency of the RGE for the isolated dijet fac- 
torization theorem, shown in Fig. [5Jb), again provides 
important constraints on its structure. Each of the func- 
tions J<[, B, and S^j includes a series of double log- 
arithms, and when these functions are convoluted over 



the kf variables at a common scale /i, these different se- 
ries have to collapse to precisely the double-logarithmic 
series of the hard function H IL . The RGE for the hard 
function H IL is a matrix equation in color space, but has 
no convolutions of kinematic variables. We expect that 
this equivalence will occur in the same manner as it does 
for the isolated Drell-Yan case. 

Key missing ingredients in providing a rigorous deriva- 
tion of Eq. ([4"9")) include i) providing a mathematically 
rigorous treatment of the separation of jets and beams in 
the factorization, and ii) determining the role of Glauber 
degrees of freedom, that in principle may couple the final- 
state jets and spoil factorization. It should be evident 
that if such a proof becomes available, it will be straight- 
forward to generalize the above discussion to the case 
where we produce N isolated jets rather than just two. 



III. THE BEAM FUNCTION 

In this section, we discuss the properties of the beam 
function in more detail. We present its definition and 
relation to the standard PDF, as well as its renormal- 
ization group evolution. We will display explicit results 
for the quark beam function at one loop (leaving a de- 
tailed derivation to a dedicated publication Ho|). The 
comparison of effects in the beam functions and PDFs 
are illustrated with plots. 

The quark, antiquark, and gluon beam functions are 
defined in SCET as 



J 



B q (Ljb + ,Lj/P-,fJ,) 

B^iob+^/p-^) 

B g {Lub+,L0/p-,fl) = -9(L0) 



Ay- 



^- e ib+ y^ 2 (p n (p-) Xn{y-^)S{co - V n ) £ X „(0) MP') 



-in 



I d V Jb+y-/2 



"{MP-) 



tr 



spin 



|xn(i/~|) S(W - V n )Xn(0)] \pn(P 



Ay_ 

47T 



„ib+y-/2 



Vn(P-) B c n ^{y-- 5(0J - V n )B%_(0) P n{P~) 



(50) 



We will briefly explain the relevant notation. (A more detailed overview of SCET and the definitions of the objects in 



Eq. (|50j) are given in Sec. II V Al ) As before, 



(l,n) and = (1, —ft) arc lightlike vectors, n 



0, n-n 



where n is a unit three- vector in the direction of the proton. The proton states \p n (P~)} have lightlike momentum 
P 11 = P~n^ /2, and the matrix elements are always implicitly averaged over the proton spin. The SCET fields for 
collinear quarks and gluons, Xn{y) and B^ ± (y), are composite fields containing Wilson lines of collinear gluons [see 
Eq. (|70j) ]. Matrix elements with these fields include so-called zero-bin subtractions [71| . which effectively divide by a 
matrix element of Wilson lines [72j ■ At lowest order in the strong coupling, the fields describe an energetic quark or 
gluon moving in the n direction with momentum p~n ,J '/2 + k^ with k <C p~ ■ The momentum operator V n picks out 
the large light-cone component p~ of all particles annihilated by Xn or &n±- Thus, when these fields annihilate the 
incoming colliding parton, the 8 function in Eq. ([50]) sets u> equal to the p~ of that parton. Therefore, x = lj/P~ is 
the fraction of the proton's light-cone momentum that is carried by the parton into the hard collision. At the time 
of the collision, this parton is propagating in an initial-state jet rather than the proton, which is encoded by the 
dependence of the beam functions on the variable b + = —k + . Here, k + = n-k is the small component of the incoming 
collincar parton's momentum. The variable t = tob + ~ — p 2 measures the parton's virtuality, where t > 0, because the 
parton is spacclikc. As we already saw in Eq. (|19p . the beam functions are convoluted with the soft function through 
b+. 
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The beam function definitions in Eq. (|50|) can be compared with those of the standard quark, antiquark, and gluon 
PDFs in SCET O, 



f q (u'/P-,n) = 6(u')(p n (P-)\x n (0) S(cj' - Vn) |xn(0)|p„(P-)) , 



fr(u'/P-,n) = 9{u')(p n (P-) 



tr sp i„ [|xn(0) <V - 7>„) X«(0)J p„(P-) 
f g (u//P-,tt) = -e(co')co'(p n (P-) Kxu{0)^-V n )B^_(0) Pn(P-) 



(51) 

The fi depend on the analogous light-cone momentum fraction £ = u)'/P~. As discussed in the beginning of Sec. Ill CI 
£ can be interpreted as the momentum fraction of the hard parton when it is taken out of the proton and due to 
perturbative corrections generally differs from x = uj/P~ appearing in the beam functions. A more common and 
equivalent definition of the PDFs is in terms of QCD fields. For example, for the quark PDF, 



f q (u'/p-,ri = 6(u>) /^e- i ^ + / 2 (p„(P-)|vi(j/+|)|w- fi ( y +|,0)^(0)|p n (p-)} 



_n\ $ 



(52) 



Equation (|51|) is essentially the Fourier transform of 
Eq. ([5"2"j) , where the SCET fields are written in momen- 
tum space with respect to u>', while the QCD fields are 
separated along the n direction between and y + h/2. 
The lightlike Wilson line Wn(y + n/2,0) is required to 
render the product of the quark fields at different space- 
time points gauge invariant, and the corresponding Wil- 
son lines in Eq. (|51| arc those hidden in the definitions 
of Xn and B£ ± . 

In the beam functions in Eq. (|50|) . the fields are in 
addition separated along the n direction, with a large 
separation y~ ^> y + corresponding to the small momen- 
tum b + -C oj. This y~ separation is formulated with a 
gauge invariant multipole expansion of fields in SCET. 
The possible gauge transformations in the effective the- 
ory are divided into global, collinear, and soft, and it is 
the coupling to soft gluons and the corresponding soft 
gauge transformations that are relevant for making the 
y~ separation gauge invariant. The collinear fields in 
Eq. ([50| are the ones that occur after making a field re- 
definition to decouple soft gluons into the soft function, 
and the resulting collinear fields no longer transform un- 
der soft gauge transformations. Hence, the SCET defini- 
tions in Eq. ([50]) are gauge invariant. 

Note that formulating equivalent definitions of the 
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FIG. 9: Tree-level diagrams for the quark beam function (a) 
and quark PDF (b). 



beam functions directly in QCD is more challenging. It 
seems to require QCD fields that are simultaneously sep- 
arated in the n and n directions, and a priori it is not 
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FIG. 10: One-loop diagrams for the quark beam function. 
Graphs (a) and (b) correspond to real gluon emission, while 
(c) and (d) are virtual corrections. Graphs (e) and (f) deter- 
mine the contribution of the gluon PDF to the quark beam 
function. 
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clear how to obtain an unambiguous gauge-invariant ex- 
pression in this case, because Wilson lines connecting 
the fields along different paths are not equivalent. For 
the beam functions, which one might think of as b + - 
dependent PDFs, this problem is solved in SCET, be- 
cause the effective theory distinguishes the large and 
small momentum components with the multipole expan- 
sion, resolving the ambiguity. 

For t = ub + ^> Aq CD , or equivalently y" <C w/Aq CD , 
the beam functions can be related to the PDFs by per- 
forming an operator product expansion in Aq CD /£ <C 1. 
This leads to the factorized form 

Z" 1 d£ / x \ 



(53) 



where j = {g, u, u,d, . . .} and 2^ is a perturbatively cal- 
culable Wilson coefficient. The physical interpretation 
of this equation was discussed in Sec. Ill CI For B g , the 
equivalent of the matching expression in Eq. (|53|) was de- 
rived in Ref. [30| for the I gg term using a moment-space 
OPE to match SCE^ onto SCET„. Ref. [IJ consid- 
ered this matching at the level of the matrix element 
defining B gi without the accompanying physical picture 
advocated here that implies that beam functions will oc- 
cur in a wide variety of interesting processes. The mix- 
ing contributions were missed in their analysis, but the 



extension of their proof to the general case is straightfor- 
ward fsoj ] . 

The coefficients Iij(t,x/£,, /i) can be determined per- 
turbatively by computing both sides of Eq. (|53[) with the 
proton states in the definitions of Bi and fj replaced by 
quark and gluon states. The tree-level diagram for the 



quark beam funct ion is shown in Fig. 9(a) and for the 
PDF in Fig. |9(b)[ They give 



B^(t,x,f,) = 6(t)S(l-x) 



from which we deduce l g r q cc (t, x/£, fi) 
In general, we have 



g»(t,| jAt )=M(*)*(i 



(54) 
S(t)5(l - x/i). 



(55) 



so the tree-level beam functions reduce to the PDFs 

B?° e {t,x,ri = 5(t)f i (x, fi ). (56) 

The one-loop coefficients for the quark beam function 
are determined from the diagrams in Fig. [TO] together 
with the corresponding diagrams for the PDFs. The 
beam functions and PDFs are renormalized using dimen- 
sional regularization with MS. From the first four dia- 
grams in Fig. [TO] we find the one-loop correction to the 
quark-quark coefficient (here z — x/£) 



2- 

5(t) 



t/n 2 



5{l-z) + - 2 

H 1 



[1-z) 



(1 - z)ln(l - z) 
1 - z 



(1 + z z ) (5(1 - z) + 9(1 - z) 1 - z 
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1 + z 2 

1 - z 
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7 

In z 



5(1 -z) 



(57) 



The last two diagrams in Fig. [TO] determine the one- loop contribution of the gluon PDF to the quark beam function, 

+ S(t) In — - 
z - 



i^ r ^ z ^ ) = ^MlLe(z)9(i- 



{ 


' 1 






J? 


[ t/v 2 \ 



[z 2 + (l-zf] +6(t) 2z(l-z) 



(58) 



At two loops, Iqq(t, z, fj,) will start to contribute as well. The plus distributions are defined as 
\e(x)g(x)} ^ = Urn 4- \0(x - e) G(x)} with G(x) = [ dx' g(x') , 
satisfying the boundary condition dx [0(x)g(x)] + = 0. In particular, 



(59) 



'dz 



[6(l-z)g(l-z)] + f(^)= ( dzg(l-z)\-f(-)-f(x)]+f(x)G(l-x) 



dt 



9(t/ii 2 )W l (t/ii 2 



_ m n +! ^ max 



n + 1 



(60) 



The infrared (IR) divergences in the diagrams in Fig. [TO] precisely cancel those in the PDF calculation as they must, 
so the matching coefficients in Eqs. (f5"T|) and ([55)) arc IR finite and independent of the IR regulator. The ultraviolet 
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(UV) divergences in the diagrams determine the one-loop RGE and anomalous dimension of the quark beam function, 
which in MS are 



d f 

fj,—B q (t,x,fJ,)= dt'^ B (t-t',fj,)B q (t',x,fJ,), J%(t, n) = -2T cusp [a s (fi)]- 



1 



Vm 2 



47T 



6C F 5(t). (61) 



Here, r cusp [a s (/x)] is the cusp anomalous dimension [74|], and the coefficient of the plus function is equal to r cusp to 
all orders in perturbation theory. The non-cusp term of the anomalous dimension is equal to that for the quark jet 
function at one-loop, and in Ref. (Hoj we prove that the anomalous dimensions for the quark beam and jet functions 
are identical to all orders in a s . As stated before, the RGE in Eq. (|6"T|) does not change x. Also, the mixing graphs in 
Figs. |l0(e)1 and [10(7)] have no UV diver gences and hence the gluon beam function does not mix into the quark beam 
functions under renormalization. Equation (1611) leads to the physical picture discussed in Sec. Ill CI The RGE has a 
solution [7(1 [lE [ill , which can be written as |67l| 



B q (t,x 7 fi)= dt' B q (t-t',x, /i ) [/ b (^',Mo,m) , U B (t, /i , A*) = e 



k b 



T(1 + Vb) 



Vb 

f4 



6(t) , (62) 



with the plus distribution defined according to Eq. (|59p . Furthermore, Kb = -ftTs(^o,/i) and r/s = T]b{ho, A*) are 



Kb(.Ho,Ij) 



da s 



4T 



cusp 



(a.) 



da' 



a s (Mo) P( a 's) 



7s(a s ) 



7?b(MO,AO 



da s 



r c usp(a! s ), (63) 



where j3(a s ) is the QCD f3 function and ~f B {a s ) is the coefficient of S(t) in j B (t,fi) in Eq. ([BTjl . 

Using the above results for the quark beam function, we can see explicitly that when we integrate over < t < i max 
to get the beam function B q (t max ,x, fj,) in Eq. (|23[) . the result contains double and single logarithms of i max //Lt , 



max ■ 

X, n) = 0(t max )f q (x, (l) + 8(t 

max J 



2tt 



C, 2 ^max 3 tmax 
F [ In — ^ In — 



A* 



+ In 



dz 

2 



0(1 -z) 



1 + z 2 
1-z 



T F [z 2 + {l-zf]f & 



(64) 



The ellipses denote ^-dependent terms that have no 
hi(imax/M 2 )- Equation (164)) shows that the natural scale 
for the beam function is \i = [Lb ~ A ma x- The loga- 
rithms of t max / fx 2 are summed by solving the beam func- 
tion's RGE in Eq. ([6"T]). From Eq. (|B"4"|) we can see how 
the matching coefficients convert the PDF running 
into the beam function running at one loop. Expanding 
Eq. (|6"Tj) to 0(a s ), the integrated beam function satisfies 

d ~ , 

(65) 

Taking the derivative of Eq. (|64|) with respect to /x, the 
first term in curly brackets proportional to f q (x,/i) re- 
produces the overall factor in Eq. (|65|) . while the terms 
in the second line precisely cancel the (i dependence of 
the tree-level term f q (x,n). Thus, even though at tree 
level B q (t,x,[i) = 5(t) f q (x, /i), the running of B q does 
not depend on x. 

To illustrate the difference between the beam func- 
tions and the PDFs, we may compare the PDFs and 
B q (t max , x, fi) for fixed i max as a function of x. For A max , 



following the discussion in Sec. Ill A[ we take 

A max = LuB+ ax = Q 2 e- 2y ^ = (xE cm ) 2 e~ 2y ^ , (66) 



1.4 r 




x 



FIG. 11: The tt-quark PDF (lower green band) and w-quark 
beam functions for y cut = 1 (middle blue band) and y cut = 2 
upper orange band. All functions are evaluated at the same 
hard scale, with the bands showing the scale variation by 
factors of two. 
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FIG. 12: The u (top row), d (middle row), and u (bottom row) beam functions at the beam scale fi% ~ f m ax for E C m = 7TeV 
and t m ax with y cut = 2 [see Eq. (|66|) ]. The panels on the left show the functions times x. The right panels show the relative 
differences compared to the respective tree-level (LO) beam functions given by the dotted lines on the left. The bands show 
the scale uncertainties as explained in the text. 



where the hard scale is taken as Q 2 = (xE cm ) 2 . By 
default we use E cm = 7 TeV and y cut = 2. For the parton 
distributions we use the NLO results of MSTW2008 [zj. 

The effect of the large logarithms contained in the 
beam function is illustrated in Fig. 1111 where we show 
the u-quark beam function xi? u (t max , x, (J-h) at NLO for 
two different values of i ma x, along with the u-quark PDF 



Hh), which is equal to the beam function at tree 
level. All solid central curves are evaluated at the com- 
mon hard scale fj, = Q = xE cm , and the surrounding 
bands correspond to varying the scale by a factor of 2. 
The lower green curve and band show the PDF (tree-level 
beam function). Including the one- loop matching correc- 
tions the beam functions are given for y cut — 1 by the 
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middle blue band and for y cut = 2 by the upper orange 
band. Clearly, with this scale choice, the beam functions 
receive large a s corrections with a dramatically increased 
scale dependence compared to the PDF, which is caused 
by the large logarithms of t max /fi 2 in Eq. (|64|) . 

To eliminate the large logarithms in the matching, 
we have to compute the beam functions at the scale 
Mb — imax; which we show in Fig. [121 Lets first con- 
sider the up-quark beam function B u shown in the top 
row. The top left panel shows xB u (t max , x, (J-b)- The 
green (light) band shows the PDF xf u (x, hh) evaluated 
at the hard scale [ah = xE cm , varying the scale by a 
factor of two (which for the top left panel is identical to 
the corresponding band in Fig. The one-loop beam 
function evaluated at fi B ~ i max is shown by the or- 
ange (medium) bands. Here, the maximum scale vari- 
ation for jjL 2 B € [imax/2, 2i max ] is not obtained at the 
edges of this region, but is closely approximated by tak- 
ing hb = {0.7, 2.0} V^max- The solid line shows the cor- 
responding central value at [ib = 1-4 y/t niax . For com- 
parison, the dotted line shows the tree-level result, i.e. 
the PDF at the scale [is- The top right panel shows 
the same curves as the top left panel, but normalized to 
the tree-level beam function. The plot shows that the 
beam function can be reliably calculated at fi = /is- Us- 
ing [i B ~ t max , the shift from the dotted to solid line 
is now of reasonable size, and the scale uncertainties are 
now similar in size to those of the PDFs. We also see 
that evaluating the PDF at the beam scale rather than 
the hard scale has a significant effect. The difference 
between evaluating the PDF at (jlh and fig is a +30% 
(—20%) correction at large (small) x, while the a s cor- 
rections from the beam function at [Ib are only ~ 10%. 
(The beam function curves increase for x > 0.4 due to 
the threshold term [ln(l — z)/(l — z)] + in l qq , but this 
region only has a small total contribution as can be seen 
on the left panel.) Since the residual [Ib dependence is 
only canceled by other contributions in the factorization 
theorem, these plots do not determine the overall size of 
the a s corrections, nor their uncertainty. These questions 
are addressed by plots of the full cross section in Scc. lVIl 

In the central and bottom panels of Fig. [T^] we show 
analogous plots for the down-quark beam function, Bd, 
and the antiup-quark, B a . While the absolute size of the 
functions for different flavors in the left panels are quite 
different, the relative corrections shown in the right panel 
are very similar. 

Fig. [12] also contains dashed lines, which show how the 
solid lines are modified if we remove the gluon contri- 
bution T qg from the one-loop beam function. The gluon 
contribution to the quark beam functions becomes no- 
ticeable at x < 0.1, increasing to about —5% at x = 0.01, 
while for the antiquark beam function it is important in 
the entire x range. This is expected, since the antiquark 
PDF is much smaller, so the gluon PDF can have a big- 
ger impact. The gluon contribution is always negative 
and partially compensates the quark matching correc- 
tion. Even at x = 0.01 there is no indication that treat- 



ing the logarithms In x in fixed-order perturbation theory 
causes any problems. The contribution from the In z term 
in Eq. (f5*7| is of similar size as other contributions and 
within the perturbative uncertainties. It has roughly the 
same size as the gluon contribution. 

IV. ISOLATED FACTORIZATION THEOREM 

In this section, we derive the isolated factorization the- 
orem in Eq. (|19p . Our analysis is based on factorization 
in SCET, which rigorously and systematically separates 
hard, soft and collinear contributions H^, [23, We 
make use of a setup with SCETi and SCETn [78j |. car- 
rying out the factorization in two stages at the scales 
Q 2 and w a ^B+ b respectively. We have an SCETj analy- 
sis to factorize initial-state jets from soft radiation. The 
initial-state jets described by beam functions in SCETi 
are then matched onto initial-state PDFs with lower off- 
shellness for the collinear particles in SCETn- In this 
section, we carry out the SCETi computation, while the 
matching onto SCETn was discussed in Sec. IIIIl Our 
analysis below uses similar tools as used in the deriva- 
tion of the factorization theorem for hemisphere invari- 
ant masses for e + e~ — > dijcts in Ref. {6^1, but differs 
significantly due to the kinematics, and the fact that we 
have initial-state rather than final-state jets and a further 
matching onto SCETn. The soft dynamics of e + e~ — > 
dijets was studied earlier in SCET in Refs. (79l[80l|. We 
start with a brief overview of the necessary SCET ingre- 
dients in Sec. IIV Al and describe the relevant kinematics 
in Sec. IIVBI We derive the factorization theorem for 
isolated pp —> XL in Sec. IIV CI including arguments to 
rule out contributions from so-called Glauber degrees of 
freedom. Finally in Sec. IIV Dl we apply the factorization 
theorem to pp — > Xt + £~ and quote final results for the 
beam thrust cross section with one-loop corrections and 
logarithmic rcsummation. 

A. SCET 

Soft-collinear effective theory is an effective field the- 
ory of QCD that describes the interactions of collincar 
and soft particles p3 - l27j . Collinear particles are charac- 
terized by having large energy and small invariant mass. 
To separate the large and small momentum components, 
it is convenient to use light-cone coordinates. We define 
two light-cone vectors 

n* = (l,n), n M = (l,-n), (67) 

with n 2 = n 2 = 0, n-n = 2, and ft is a unit three- vector. 
Any four-momentum p can then be decomposed as 

P M = n-p— + n-p — +p^ x - (68) 

Choosing n close to the direction of a collinear particle, 
its momentum p scales as (n-p, n-p,p n ±) ~ h-p (A 2 , 1, A), 
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with A <C 1 a small parameter. For example, for a jet 
of collinear particles in the ft direction with total mo- 
mentum px , n ■ px — 2Ex corresponds to the large 
energy of the jet, while n-px — Px/Ex -C Ex, so 
\^ Px /E x «l. 

To construct the fields of the effective theory, the mo- 
mentum is written as 



n-p — 

y 2 



(69) 



where n ■ p ~ Q and p n j_ ~ XQ are the large momentum 
components, where Q is the scale of the hard interaction, 
while k ~ X 2 Q is a small residual momentum. The effec- 
tive theory expansion is in powers of the small parameter 
A. 

The SCET fields for n-collinear quarks and gluons, 
£n,p(%) and A n ^{x), are labeled by the collinear direction 
n and their large momentum p. They are written in po- 
sition space with respect to the residual momentum and 
in momentum space with respect to the large momentum 
components. Frequently, we will only keep the label n de- 
noting the collinear direction, while the momentum labels 
are summed over and suppressed. Derivatives acting on 
the fields pick out the residual momentum dependence, 
id" ~ k ~ X 2 Q. The large label momentum is obtained 
from the momentum operator V^, e.g. Vfi £ n> p = p" £, n ,p- 
If there are several fields, V n returns the sum of the la- 
bel momenta of all n-collinear fields. For convenience, 
we define V n = n-V n , which picks out the large minus 
component. 

Collinear operators are constructed out of products of 
fields and Wilson lines that are invariant under collinear 
gauge transformations [25|, |2(|. The smallest building 
blocks are collinearly gauge-invariant quark and gluon 
fields, defined as 



Xn,M)= [6(u-V n )Wl(x){; n (x) 

K, u ±(%) = ^[s(w + V n )Wkx)U% x W n (x) 
where 

is the collinear covariant derivative and 



W n (x) 



perms 



(70) 



(71) 



(72) 



The label operators in Eqs. ([70]) and (J72J) only act inside 
the square brackets. Here, W n (x) is a Wilson line of 
n-collincar gluons in label momentum space. It sums 
up arbitrary emissions of n-collinear gluons from an n- 
collinear quark or gluon, which arc 0(1) in the power 
counting. Since W n (x) is localized with respect to the 
residual position x, we can treat Xn,u( x ) an d B^ u (x) as 
local quark and gluon fields. The label momentum to is 
treated as a continuous variable, which is why we use a 6- 
function operator in Eq. (fT0"|) . It is set equal to the sum of 



the minus label momenta of all fields that the 5 function 
acts on, including those in the Wilson lines, while the 
label momenta of the individual fields arc summed over. 

In general, the effective theory can contain several 
collinear sectors, each containing collinear fields along a 
different collinear direction. To have a well-defined power 
expansion in this case, the different collinear directions 
m have to be well separated [73| . 



» A 



for 



i ^ 3 , 



(73) 



which is simply the requirement that different collinear 
sectors are distinct and do not overlap. For pp — > X£ + £~ , 
we need two collinear sectors, n a and rib, along the di- 
rections of the two beams. We use a bar to denote the 
conjugate lightlike vector, so n, ■ = 2. As the beams 
are back-to-back, we have n a ~ rift, so n a -ni, ~ 2 and 
Eq. (jT3"|) is easily satisfied. 

Particles that exchange large momentum of 0(Q) be- 
tween collinear particles moving in different directions 
have to be off shell by an amount of C(n, • rijQ 2 ). These 
modes can be integrated out of the theory at the hard 
scale Q by matching full QCD onto SCET, which yields 
the hard function. The effective theory below the scale 
Q then splits into several distinct collinear sectors, where 
particles in the same collinear sector can still interact 
with each other, while at leading order in the power 
counting particles from different collinear sectors can 
only interact by the exchange of soft particles. This 
means that before and after the hard interaction takes 
place, the jets described by the different collinear sectors 
evolve independently from each other with only soft but 
no hard interactions between them. 

The soft degrees of freedom, responsible for the radi- 
ation between collinear jets, arc described in the effec- 
tive theory by soft 10 quark and gluon fields, q s (x) and 
A s (x), which only have residual soft momentum depen- 
dence i9 M ~ X 2 Q. They couple to the collinear sectors 
via the soft covariant derivative 



wt 



id" + gA^ 



(74) 



acting on the collinear fields. At leading order in A, n- 
collinear particles only couple to the n-A s component of 
soft gluons, so the leading-order n-collinear Lagrangian 
only depends on n-D s . For n-collinear quarks [25l |26| 

C n = Um-D s + gn-An + ip n± W n ^- W n i%A fen • 

(75) 

The leading-order n-collinear Lagrangian for gluons is 
given in Ref . [27| ■ 



10 In some situations it is necessary to distinguish two types of soft 
sectors, referred to as soft and ultrasoft in the SCET literature. 
In this paper we only need what are usually called ultrasoft par- 
ticles, so we will simply refer to these as soft. 
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The coupling of soft gluons to collincar particles can 
be removed at leading order by defining new collinear 
fields H3 

X { nl( X )= Y n( X )XnA^), (76) 

where Y n (x) and y n {x) are soft Wilson lines in the fun- 
damental and adjoint representations, 



Y n {x) = Pexp 



ig / ds n-A s (x 

J — oo 



T c y c n d {x) =Y n {x)T d Y^x) 



(77) 



The symbol P in Eq. (|77jl denotes the path ordering of 
the color generators along the integration path. The in- 
tegral limits in Eq. (|77[) with the reference point at — oo 
are the natural choice for incoming particles (8i| . The 
final results are always independent of the choice of ref- 
erence point, and with the above choice the interpolating 
fields for the incoming proton states do not introduce 
additional Wilson lines |82| . 

After the field redefinition in Eq. (|76|) . the leading- 
order SCET Lagrangian separates into the sum of inde- 
pendent rii-collinear and soft Lagrangians, 



£scet — ^ C$ + £s 



(78) 



with no interactions between any of the collinear and soft 
sectors. The ellipses denote terms that are subleading in 
the power counting. This decoupling is what will allow 
us to factorize the cross section into separate beam and 
soft functions. The field redefinition in Eq. (|76|) intro- 
duces soft Wilson lines in the operators, which because 
of Eq. ([78]) can be factored out of the matrix element and 
will make up the soft function. 



the collision, and are also responsible for initiating the 
hard interaction. We define the momenta of the spacelike 
off-shell partons that go into the hard interaction as 



2 



It',. 



p b — x b & cm — b — 0,, 



b_L ' 



(80) 



where x a and Xb are the light-cone momentum fractions 
at which the beam functions will be evaluated. The 
power-counting parameters for the collinear sectors are 



A 



b+ 



x b E c 



(81) 



where the relevant momenta are those of the off-shell par- 
tons in Eq. (|80|) . because these are the momenta carried 
by the n a - and 7i&-collmear fields. 

We write the momentum of the incoming partons that 
are taken out of the proton as 



^cm^ +0(A QCD ). 



-O(Aqcd), (82) 



which defines the light-cone momentum fractions £ aj b at 
which the PDFs are evaluated. The typical _L-momcnta 
of partons in the proton are O(Aqcd), while the small 
plus components arc 0(Aq CD /i? cm ). These momenta 
are much smaller than any soft or residual momenta in 
SCETi and arc expanded, which precisely corresponds to 
the OPE for the beam functions in Eq. (|34|) when match- 
ing them onto SCETn. 

The momentum of the final-state remnant of the pro- 
ton is thus given by 



(1 - ta)E c 



2 ' 



(83) 



while the remnant of the initial-state jet radiated into 
the final state by the beam function has momentum 



B. Kinematics 

Before deriving the factorization theorem, we discuss 
the relevant kinematics, as illustrated in Fig. [T3J As al- 
ready mentioned, we introduce a separate set of collinear 
fields for each of the beams, with the light-cone vectors 
n a and n b aligned with the beam directions. To derive 
the factorization theorem we work in the center-of-mass 
frame of the hadronic collision, so the momenta of the 
incoming protons are (neglecting the proton mass) 



Pff = E r 



pi' 



E r 



(79) 



with n a = —fib- In particular, rib = n a and n a -rib = 2. 
We will mostly keep the dependence on the two beam 
directions explicit, but one should keep in mind that n a 
and rib are related. 

The collincar fields in the n a and rib directions describe 
the interactions within each of the beams before and after 



b» = (& - x a )E cm ^+ b +^+b^, (84) 

and similarly for the rib direction. The total n a -collinear 
momentum in the final state is the sum of Eqs. (|83p and 
(|84|) . or equivalently, the difference between the proton 
momentum and Eq. (|80|) . 

K+rZ = n-V^a = (l-x a )E cm -f+bt -f+b^ . (85) 

In addition to the collinear momenta, we define k 1 ^ as 
the total four-momentum of the soft radiation in the final 
state. Hence, the total hadronic momentum in the final 
state is given by 



P' 1 X = (P^-P^ + (Pb"-p"b) + k 



S 1 



(86) 



and we can write total momentum conservation 
Pb = Px + ^ as 



Pt+P^q^ + k^ 



(87) 
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FIG. 13: Definition of the different collinear momenta related to the incoming beams. The soft radiation is not shown. 



where is the total leptonic momentum. 

The collinear and soft momenta, b%, k£ are not ex- 
perimentally measurable quantities. Instead, the experi- 
ments can only measure hadronic quantities, such as the 
hemisphere momenta = n a ■ B a and B b — n b - B b in- 
troduced in Sec. Ill Al Splitting the total soft momentum 
into its contributions from each hemisphere, k% = k£ + k b 
as shown in Fig. [3l we then have 



B£ = K + r£ + k£ 



= + < + (88) 



and defining fc+ = n a -k a , k b = n^-kb, we get 

B+ = n a -B a = b+ + k+ , B+ = n b -B b = b+ + k+ . 

(89) 

In particular, the remnant momenta b do not con- 
tribute to B^ b . A physical argument for this was dis- 
cussed in Sec. Ill Al 

Next, we decompose the total leptonic momentum as 



to the transverse components qr such that we can freely 
integrate over them. Therefore, we have to expand the 
kinematics in the limit qr = 0. This expansion is justified 
because from Eq. (18TB we have 



~Px± 



-Kx-K 



b± 



k^ 
K s± 



XQ ■ (94) 



A parametrically large q^_ ~ Q would require a separate 
jet at large px ~ Q to balance the transverse momentum, 
which is not allowed in our setup. After integrating over 
the lepton phase space this expansion incurs power cor- 
rections of order q\ ~ \ 2 Q 2 . The kinematics of the hard 
matrix element in the factorization theorem is then given 
by the tree-level partonic kinematics, with the partonic 
momentum conservation 



•Ea-B c 



which implies 



x b E c 



-K 4 

q = q f +q T 



(95) 



fi — a 



I 1 



q-q -f + q 2 +q± ' 



(90) 



where q^_ contains the two components of g M transverse 
to the beam direction. Taking the z-axis along the ft a 
beam direction, we have 



q ± = q" T q z , 



<# = (0,fr,0) 



(91) 



where qr = {q x ,q y ) is a two- vector in the transverse x- 
y-plane. The total leptonic invariant mass and rapidity 
are 



x a E, 
x b E, 



cm 


= q = 


Vq e , 


'cm 


= q + = 


r~2 -y 


q 2 


= q + q~ 




Y 


= Ilni 

2 q 


— = — In — 

f 2 x b 



(96) 



Equations (|94|) and (|95|) imply that parametrically the 
leptons are back-to-back in the transverse plane. Since 
q + and q~ can differ substantially, the leptons do not 
need to be back-to-back in three dimensions. 



2 + - , 2 + - -2 

q =q q +q± = q q -qr, 



with 



F=lln^-, (92) 



d 4 q = - dq+dq- d 2 q T = - dq 2 dY d 2 q T 



(93) 



As we will see in the next subsection, the derivation 
of the factorization theorem requires us to be insensitive 



C. Derivation of Isolated Factorization Theorem 

We now proceed to derive the isolated factorization 
theorem for generic processes pp — > XL, where the 
hadronic final state X has a restriction on the hemi- 
sphere momenta B^ b . The derivation is carried out using 
SCET without Glauber degrees of freedom. The proof 
that Glauber effects are not required is given at the end 
of this subsection. 



20 



1. Cross Section in QCD 

We will generically refer to properties of L as "lep- 
tonic" , even though L can contain any non-strongly in- 
teracting particles. We only consider processes where 
the hard interaction couples the strong and electroweak 
sectors through one two-particle QCD current. (This 
includes for example Drcll-Yan or Higgs production 
through gluon fusion with the Higgs decaying non- 
hadronically, but does not include electroweak Higgs pro- 
duction via vector-boson fusion.) Then, at leading or- 
der in the electroweak interactions, we can factorize the 
full-theory matrix element into its leptonic and hadronic 
parts 



M (pp -> XL) = J2 L J ( x | j \pp) 



(97) 



The sum runs over all relevant color-singlet two-particle 
QCD currents J, and Lj contains the corresponding elec- 
troweak matrix element, including the electroweak prop- 
agator coupling to J. For example, for Drcll-Yan with 
L = £ + £~ , the relevant currents are 



J Vf = 4f^<lf> 



J Af 



9/7 75<7/ , 



(98) 



so in this case the sum over J in Eq. (|97[) includes the 
sums over the two Dirac structures, the vector index /z, 
and the quark flavor / = {it, <i, . . .}. The corresponding 
Lyj: and L^ are given below in Eq. (|137p . 

The cross section for some hadronic observable O in 
the center-of-mass frame of the collision, averaged over 
proton spins, is 



da 



1 



d 2 q T 



dq 2 dYdO 2E 2 m J 2(2tt) 4 



spins X 
\4r4 



x (2nY8\P a + P b -q~p x ). 



(99) 



Here, P a , b are the incoming proton momenta, px and 
PL are the total hadronic and leptonic momenta, d$i 
denotes the leptonic phase space, and the phase-space 



integrations for the hadronic final states are included in 
the sum over X. The last S function is overall momen- 
tum conservation. The function fo(X) inside the second 
5 function returns the value of the hadronic observable 
O for a given hadronic state X, so the S function picks 
out all final states that contribute to a certain value of 
O. The S 4 (q — pl) under the leptonic phase-space inte- 
gral defines the measured q as the total leptonic momen- 
tum. Expanding this 5 function for qx — 0, the leptonic 
part does not depend on qr at leading order, and using 
Eq. (|97|) . we can rewrite Eq. (|99|) as 



da 



d<7 2 dYdO 2£ c 2 m 



^^L^(g 2 ,F)^(g 2 ,y,0). 

(100) 



The leptonic tensor is defined as 
Lj. r (q 2 ,Y)=Jd<S> L L\Lj, (2n) 4 S 4 (c 



rib 



PL 



2 * 2 

^ (101) 

where q ± = \fq^e^ Y . The hadronic tensor contains the 
square of the hadronic matrix element 



Wjj,(q 2 ,Y,0) 
d 2 q T 



■Jpp\jH0)\X)(X\j'(0)\pp) (102) 

x (27r) 4 <5 4 (P a + P b - q -px) S[0 - fo(X)} , 

where as in Sec. lIIII we keep the average over proton spins 
implicit in the matrix element. Since Wjji is integrated 
over <fr, it can only depend on q 2 and Y, as well as the 
hadronic observable O. 

We are interested in the hadronic observables £>+ = 
n a ■ B a and — n b ■ B b . The hemisphere hadronic 
momenta B^ b (X) can be obtained from the states \X) 
using the hemisphere momentum operators p 1 ^ b 

ft\X) = B£{X)\X) , %\X) = Bg(X)\X) . (103) 

A ficld-thcorctic definition of p^ b in terms of the energy- 
momentum tensor of the field theory was given in 
Ref. [13 ■ The hadronic tensor for O = {B+ , S+} is 



Wjj,(q 2 ,Y,B+,B+) = J r ^xe-' lq - x Y.(M J ^ x )\ X )( X \ J '^)\PP) 5 [ B i - n a - B a (X)] 5[B+ - n b - B b (X)] 



dx + dx 
(4tt) 2 



Ji(x-'f +x + '-f)8(B+ - n a -p a )S(B+ - n b -p b )J'(0) 



PP. 
(104) 



In the first line we used momentum conservation to shift the position of , and in the second line we performed the 
integral over gV, which sets xt to zero. We also used Eq. (|103|) to eliminate the explicit dependence on X, allowing 
us to carry out the sum over all states X. The restriction on the states X is now implicit through the operator 5 
functions inside the matrix clement. 
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2. Matching QCD onto SCET 



In the next step, we match the QCD currents J onto SCET currents by integrating out fluctuations at the hard 
scale Q. At leading order in the power counting, the matching takes the form 



J(x) 



duii do;2 e 



-i(bi+b 2 )-x 



E c t^ o a J(biM- x) + c^ihM) o gg ^(b u b 2 ; x) 



(105) 



where a, j3 are spinor indices, /i, ;/ are vector indices, and the sum over q runs over all quark flavors {u, d, . . .}. The 
Wilson coefficients and operators depend on the large label momenta 



(106) 



They will eventually be set to either q~n£/2 or q + n^/2 by momentum conservation, but at this point are unspecified, 
and the sums and integrals over rt\ , ri2 and wi , W2 in Eq. (1 10511 run over all sets of distinct collinear directions and 
large label momenta. On the right-hand side of Eq. (|105p . the full x dependence of the current is separated into the 
x dependence appearing in the overall phase factor with large label momenta and the residual x dependence of the 
SCET operators. 

The SCET operators O gq ^(x) and 0^ g (x) arc constructed out of the collinear fields in Eq. ([70|) . At leading order 
in the power counting they contain one field for each collinear direction. Since the QCD currents are color singlets, 
the leading operators that can contribute are 



0*nbi,b 2 ;x) 



(107) 



where j and c are color indices in the fundamental and adjoint representations. We included appropriate minus signs 
on the labels, such that we always have u>i t 2 > for incoming particles. Here, x = Xq is a quark field of flavor q, 
which for simplicity we keep implicit in our notation. Note that the entire spin and flavor structure of the current J is 
hidden in the label J on the matching coefficients in Eq. (|105|) . The gluon operator is symmetric under interchanging 
both fx v and b\ o 62, so its matching coefficient must have the same symmetry, 



We define the conjugate quark operator and matching coefficient with the usual factors of 7 , i.e., 



Ojf a (6i,6 2 ) = X^MxZi,-^) . ^(61,62) = [7 U ^-(6 1 ,6 2 )7 



(108) 



(109) 



The matching coefficients are obtained by computing the renormalized matrix elements (0|...|gg) and (0|...|<?<?) on 
both sides of Eq. (|105|) and comparing the results. In pure dimensional regularization for UV and IR divergences all 
loop graphs in SCET are scaleless and vanish, which means the UV and IR divergences in the bare matrix elements 
precisely cancel each other. The renormalized matrix elements of the right-hand side of Eq. (j!05p are then given by 
their tree-level expressions plus pure 1/e IR divergences, which cancel against those of the full-theory matrix elements 
(0|J|gg) and (0| J\gg) of the left-hand side. Hence, the matching coefficients in MS are given in terms of the IR-finitc 
parts of the renormalized full-theory matrix elements computed in pure dimensional regularization. 



3. Soft-Collmear Factorization 
The field redefinitions in Eq. (|76|) introduce soft Wilson lines into the operators in Eq. (|107l) , 

Ot (x) = xT-^ (*) T [Yl (x) Y n2 (x)] j \^ 2 (x) , 

0%(x) = >/^b£^ Ui± (*)t^^ . (110) 

The time ordering is required to ensure the proper ordering of the soft gluon fields inside the Wilson lines. It 
only affects the ordering of the field operators, while the ordering of the color generators is still determined by the 
(anti)path ordering of the Wilson lines. In the remainder, we use these redefined fields and drop the (0) superscript 
for convenience. 
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Since the momentum operator is linear in the Lagrangian, Eq. (|78|) allows us to write the hemisphere momentum 
operators as the sum of independent operators acting in the separate collinear and soft sectors, 



Pa Pa,n a T" Pa,n b "T Pa,. 



Pb = Pb,n a + Pb,n b + Pb,s 



(111) 



The n a (n b ) collinear sector cannot contribute momentum in the n b (n a ) hemisphere. Thus, p a ,n b = Pb,n a = 0, while 
Pa,n a = Pn a and p*b,n b = Pn b reduce to the total momentum operators for each of the collinear sectors. For the soft 
sector, the distinction between the two hemisphere operators is important. We can now write 



- h t - k t) S ( b a - n a-Pn a )S(k+ - n a -p a ^) , 



5{B+-n a -p a ) = Jdb+dk+S(B+ 
6(B+ - n b -p b ) = Jdb+ dk+ 5{B+ - b+ - k+) S(b+ - n b -p„ b )5(k+ - n b -p b . s ) , 



(112) 



Using Eq. (|105p in the hadronic tensor in Eq. (|104p . the forward matrix element of the product of currents turns into 
the forward matrix element of the product of the operators in Eq. (jllOp . Since the Lagrangian in Eq. (|78p contains 
no interactions between the collinear and soft sectors after the field redefinition, we can use Eq. (I112p to factorize the 
resulting matrix element into a product of independent n a -collincar, nf,-collinear, and soft matrix elements. 
We first look at the contribution from O qq . The x integral of the forward matrix element of O qq becomes 



dx^dx' 



(4tt) 2 
' dx+ dx 



„ i(a+ X - +q - X + )/2 e i(b 1+ ~b 2 )- X {Mnb 



0\l a {x) S(B+ - na-Pa) 5(B+ - n b -pb) O a q { (0) 
e - Kq+x - +q - x+)/ 2 J d6 + db + dk + dk + §{B + _ b +_ k + } S{B + _ b +_ k + } 



Pn a Pn b 



(47T) 

du a du b e ^^ + +^bx )/2 J fi( U2 _ U( j _ _ ^ (j^/ _ 



x e(u a )( Pna x P n k (x) S(b+ - n a -p na )8{u a - V na )xC (°) 



/3' k' i 



Pi 

Pn t 



X (0 



T[Yl(x)Y nb (x)] kj S(k+ -n a -p a , s )6(k+ - n b -p b ^)T[Yl(0)Y na (0)Y O) + (a o b) 



■ i , i 



(113) 



Here, \p na ) and \p Ub ) are the proton states with momenta P^ b = E cla n^ L b /2 as in Eq. (|79p . The two terms in brackets 
in Eq. ()113p arise from the different ways of matching up the fields with the external proton states. The restriction to 
have positive labels w requires the fields in O qq to be matched with the incoming proton states and the fields in O qq 
with the outgoing proton states. In principle, there are two more ways to match the fields and external states, yielding 
matrix elements with the structure (p|xx|p) and (p|xx|p), which vanish due to quark flavor number conservation in 
QCD. For the same reason, in the full product (J^ O qq )(J2 q ' Oq'q') onr y t ne flavor-diagonal term with q = q' survives. 

We abbreviate the collinear and soft matrix elements in the last three lines of Eq. (|113p as M UJa (x^), M UJb (x + ) 1 
M s (x + , x~). The collinear matrix elements only depend on one light-cone coordinate because the label momenta uj a ^ b 
are defined to be continuous. We could have also started with discrete label momenta, Cj a ,b, and then convert to 
continuous labels by absorbing the residual k~ dependence as follows: 



c+ / 2 M £o (*+,*-) 



E 



dk; 



C+ / 2 A/ Wa (x-), (114) 



and analogously for M Ub (x + ). In the second step we used that by reparametrization invariance the Fourier-transformed 
matrix element can only depend on the linear combination uj a + k~ = ui a . 

As an aside, note that in the well-studied case where the collinear matrix elements are between vacuum states, 
giving rise to jet functions, the distinction between discrete and continuous labels is not as relevant. In that case, the 
SCET Feynman rules imply that the collinear matrix elements do not depend on the residual k~ (and k±) components, 
and therefore the label momenta can be treated in cither way. In our case, momentum conservation with the external 
state forces the collinear matrix elements to depend on k~ . Therefore, the only way to eliminate the residual k~ 
dependence is to absorb it into continuous to labels. One can easily see this already at tree level. Replacing the proton 
states by quark states with momentum p = p + p ri we get 



' dx + dx 
(4tt) 2 



\(k~x + +k + x~)/2 



(q(p)\xn(x + , x )S Q 7p n x„(0)\q(p)) =uu5 SlS -8(k - p r ) 5(k + - p+) . 



(115) 
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and the label and residual minus momenta are combined using 8^ p- 5(k~ — p~) = 5(u> — p~)- Our continuous u> is 
physical and corresponds to the momentum fraction of the quark in the proton. 

Returning to our discussion, to perform the x integral in Eq. (| 1 13[) . we take the residual Fourier transforms of the 
matrix elements, 



M Ua (x~) = 
M s (x~,x + ) = 



AU+ . 

e ^*-/2 Moja{k + h Mcjb{x+) 



2tt 
dfc+ dfc 



dk_ 



jk-x + /2 Mub{k -^ 



fL e Kkt* + +K* + )/* Ms (kt,kj). 



(116) 



Just as x^ , the residual momenta k^ and kf here are all defined with respect to the common n — n a . The x integral 
in Eq. (|1 X3[) now becomes 



dx + da; 



J{u a - q -)x + /2 e i(u„-q+)x 72^^-) M Ub ( X + ) M S (X+,X-) 



dk+ dfc" dfc+dfc: 



M Wo (fc + ) M Ub (k-)M.(k+,kJ) 5(oj a -q-+k- + kj) S(cj b - q+ + k+ + fc+) 



2tt 2tt (2tt) 2 
6(u a - q-)S(u b - q + ) M Ua (0) AL b (0) M.(0) 



(117) 



In the last step we expanded q ,± — fc* — kf = q ± [l + 0(A 2 )]. The remaining residual integrations are then simply the 
Fourier transforms of the matrix elements at x = 0. 

Note that without the integration over qr in the hadronic tensor Eq. (|104p . the currents would depend on x±> 
which would require us to include perpendicular components 6 a ,&± in the label momenta, and the soft matrix element 
would depend on x±, too. (The residual fcj_ dependence in the collinear matrix elements can again be absorbed into 
continuous b a>b ±.) The corresponding x± integration in Eq. (| 1 1 7[) would yield an additional 5 function 8 2 (b a ± + b b ± + 

qr — k s ±). Integrating over qr effectively eliminates this 8 function, which would otherwise force us to introduce an 
explicit dependence on b a , b ± in the beam functions. If one considers the qr spectrum of the dileptons for q\ <C q 2 , 
our analysis here provides a starting point but requires further study. One cannot just use pr in place of B^ b with 
our arguments to impose an analogous restriction on the final state, because at 0{a 2 s ) one can have two jets at high 
Pt that still have small total pr- 

The n a -collinear matrix element now reduces to the quark beam functions defined in Eq. (|50D , 



M Ua (0) = 0{u, a ) {p na xi k a (0) S(b+ - n a -p na ) 8(uo a - V na ) xl* (0) 



/3'fe' 



if 


8 k ' k 


4 






S k'k 


4 






S k'k 


4 





Y 6 lle{uj a )(p na \%n a (0) 8(b+ - n a -p na )S(u a - V na ) ^ X n a (0) \ Pna ) 



4-7T 



e -if>i a y /2 e ipt a V /2 



Xn.(0)< 



2 



(118) 



We abbreviated p+ a = n a -p na , and in the last step we used e lj3 ™ y ^ 2 Xn(0)e lp ™ v I 2 = Xn(y n/2) andp+|p n ) = 0. 
Similarly, for the antiquark beam function we have 



M„M = 0(u b )( Pnb j<%(x) 8(b+ - n b -p b ) 8(u> b ~ P„ b )Xn/(0) 



Pr, 



H a ' 8 jj ' 

4 Nr. 



uj b Bg(to b b+,uj b /P b ) 



(119) 



Since the collinear matrix elements are color diagonal, the soft matrix element reduces to an overall color-singlet 
trace, which defines the qq incoming hemisphere soft function, 

SZ mi (k+,k+) = i- tr (0\T[Yl (0) Y nb (0)] 8(k+ - n a -p a , s ) 8(k+ - n b -p b , s )T[Yl (0) Y Ua (0)] |0> . (120) 

The trace is over color and the factor of 1/N C is included by convention, such that at tree level we have 
SihcraT (k£ , k b ) — S(k^) S(k b ). The soft matrix element in the second term of Eq. (|113|) with a O b interchanged 
is equal to the above one due to charge conjugation invariancc of QCD. Under charge conjugation, the Wilson lines 
transform as C~ l Y£>C = TfY^p'*]. The explicit time ordering is required because the fields in Y n are time-ordered 



30 



by default, and charge conjugation only changes the ordering of the color generators but not of the field operators. 
For us this is not relevant, because the ordering of the fields is determined by the overall (anti-)time ordering in the 
matrix element. Thus, for the soft matrix element with a «-» b interchanged, we find 



tr(0\T[Y% b Y na ]8(k+ - n a -p a , s )8(k+ - n b -p b , s ) T[Y^Y nb ] |0) 

£ to(0\T[Y%Y£]6(k+ n a -p a , s ) 8(k+ n 6 -p 6>s ) T^F^] |0) = 5^ mi (fc+, k+) 



(121) 



where the transpose refers to the color indices. In the last step we used tr[A T B T C T D T ] = tr[BADC] and the fact 
that the fields in Yj[ and Y na are spacclikc separated and thus commute. Under parity, we have P _1 y„ a P = Y nb and 
P n a -p a)S P = n b -p bt3 . Therefore, CP invariance implies that S?£ ■ is symmetric in its arguments, 

S$Li(k+,k+) C = P ±- tr (0\T[YlY nb ]S(k+ - n b -p b , s ) 8(k+ n a -p a , s ) T[Y^ h Y na ] |0) = S^ mi (A+ fc+) . (122) 

Having worked out the different terms in Eq. (| 113[) . we are ready to include the remaining pieces from Eqs. (|105p 
and (|104p . The qq contribution to the hadronic tensor becomes 



Wj rq ,(q\Y,Bt,B+) 

du a duj b 8 (u) a - q~ ) 8(w b - q + ) ^ 

ni ,ri2 ,ti'^ ,n' 2 



dwi dw 2 doj[ dio' 2 c%(b x , l 2 ) c a / Q _ (y x , y 2 ) 



j'<a y 



x \ S n 2 n a 5(UJ 2 - Wo) 5n' 2 n a H^' 2 ~ W a ) <5«m b 6(u)l - UJ b ) 5 n ' lTlh 8(u[ - LxJ b ) (123) 

i- J dfc+ dk+ q 2 B q [x a E cm (Bt - k+), x a ] B q [x b E cm (B+ - k+), x b ] S?,f cmi (fc + , k+) + (a o b) J 



= Hjj, qcj (b a ,b b ) J dfc+dfc+ q'B q [x a E cm (B+ - k+),x a ] B q [x b E cm (B+ - k+), x b ] 5f h 9 cmi (fc+, fc+) + (g O g) . 
All label sums and integrations from Eq. (| 105[) eliminate the label J's from Eq. (|113|1 . In the second step we defined 



rp a_ rp o_ 



E c 



E c 



/ q z e 
E cm 



x b 



u b 



9 — y 
g z e 

Er.m 



as in Eq. (|96|) . and introduced the hard functions 



Hjj> qq -(b a , b b ) — — - tr spins Cj qq (b b , b a ) Cj> qq (b b , b a ) 



Hjj' qq (b a ,b b ) — Hjj' q q(b b ,b a ) 



(124) 



(125) 



Equation (|123|) is the final factorized result for the O qq contribution to the hadronic tensor. 
Repeating the same steps for O gg , we obtain for the forward matrix element 



dx + dx 
(4tt) 2 



-i{q + x-+q-x + )/2 „i(&i+fc 2 )x 



Pn a Pn„ 



Ol? (x) 6(B+ n a -p a , s ) 8{B+ n b -p b , s ) off (0) 



Pn a Pn b 



du a d^ h 8{Lo a - q~) 8(cj b - q+) J db+ db+ dk+ dk+ S(B+ - b+ - k+) S(B+ - b+ - k+) 
8(oj 2 - u b ) + (a b)] [8 n ' i1la 8(u[ - us a ) 8 n '^ nb 8(uj 2 - uj b ) + (a «-» 
x w a ^(cj ) (p„ a Bf x (0) <J(6+ - n a -p„J <5K - P„ a ) S^JO) 



x Wb ff( W6 ) ( Pnb B^ ± (0) 5(6+ - %6 - ^nJ^A(°) 



x (0 



T [yJL (°) ^ (0)] Cd S(k+ - n a -p a , s ) 8(k+ - n b -p 6>s ) T (0) y„„ (0)] 



(126) 



where we already performed the integral over x. The four terms in the third line correspond to the four different 
ways to match up the gluon fields with the incoming proton states. The collinear matrix elements reduce to the gluon 
beam function defined in Eq. (|5"0|) , 



u a 6(oj a ) (p na B£ c _l(0) 8(b+ - n a -p na )S(oj a - V„ a ) c x (0) 



Pn a ) = 



8 CC 
2 W - 1 



^fl ff (w„6+,w /P-). (127) 
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Including the color traces from the beam functions, the soft matrix element defines the gluonic incoming hemisphere 
soft function, 

SZmi(ka,kt) = J^r-J (o|tr color {T[yt a (0)y„ b (0)] 6(k+ - n a - Pa , s )S(k+ - n b -p b , s ) T [y\ b (0) y na (0)] } |o) , (128) 

where the normalization is again convention. Putting everything together, the gluon contribution to the hadronic 
tensor becomes 

Wjj, gg (q 2 , Y, B+,B+) = Hjj, gg (b a , b b ) J dfc+dfc+ q 2 B g [x a E cm (B+ - k+),x a }B g [x b E cm (B+ - k+),x b ] 5f h » mi (A:+, k+), 

(129) 



with the hard function 



H. 



1 1 



JJigg(ba,b b ) = j^f- - ^ (g_L W ' 9Xvu> + 9± »u> 9X u»' ) Cj^* (b a , b b ) C% v gg (b a , b b ) . (130) 



Here we have used the symmetry of the Wilson coefficients in Eq. (|108p to simplify the four terms that arise from 
interchanging a -f-s- b in Eq. (|126j) . 

To obtain the full result for the hadronic tensor all we have to do now is to add up the contributions from the 
different quark flavors and the gluon, 

Wjj, (q 2 ,Y, B+,B+) = ]T Wjj, qq -(q 2 , Y, B+,B+) + Wjj, gg (q 2 , Y, B+,B+) . (131) 

Inserting this back into Eq. (| 100[) . the final result for the factorized cross section becomes 
da 



dq 2 dYdB+dB u 



J2H l3 (q 2 ,Y) f dkidk+q 2 B l [x a EUBt-kt),x a ]B ] [x b E cm {B+-k+),x b ]Sll^ (132) 

6 -- J 



with x atb E cm = 
the hard function 



q 2 e ±Y as in Eqs. ([96]) and (fT24l) and 



Hi 



j, j' 



Hjj'i 3 (xaE cm '-^-,x b E cm ^-^ . (133) 



The sum in Eq. (|132|) runs over parton species ij = 
{gg,uu,uu, dd, . . .}, where Bi is the beam function for 
parton i in beam a and Bj for parton j in beam b. Equa- 
tion (jl32[) is the final factorization theorem for the iso- 
lated pp — > XL and pp —> XL processes. In Sec. IIVDI 
below we will apply it to the case of Drell-Yan, which will 
yield Eq. ([15)1 . 

The beam functions in Eq. ()132p are universal and 
take into account collinear radiation for isolated processes 
with x away from one. Since the soft function only de- 
pends on the color representation, but not on the specific 
quark flavor, there are only two independent soft func- 
tions SX, mi and S^emi- I n the sum over ij i n Eq. (|132D . 
there are no mixed terms with ij corresponding to beam 
functions of two different quark flavors. Likewise, there 
are no mixed terms with quark and gluon beam func- 
tions. For example, a graph like Fig. 6(e) is part of the 
ij = qq term in the sum. Thus, cross terms between 
quark and gluon PDFs only appear via the contributions 



of different PDFs to a given beam function, as shown in 
Eq. ([341). 

The only process dependence in Eq. (|132p arises 
through the hard functions Hij(q 2 ,Y), and one can 
study any desired leptonic observables by inserting the 
appropriate projections in the leptonic phase-space in- 
tegrations inside Ljj'(q 2 ,Y). Since the hard func- 
tion Hjj>ij corresponds to the partonic matrix element 
(ij\ J' i '|0)(0| J'\ij) and Ljji is given by the square of the 
relevant electroweak matrix elements L'jLji, Hij(q 2 ,Y) 
can be determined from calculations of the partonic cross 
section ij — > L. Furthermore, Hij(q 2 ,Y) is identical to 
the hard function in threshold factorization theorems and 
hence in many cases is known from existing computa- 
tions. 



4- Cancellation of Glauber Gluons 

In the above derivation we have implicitly assumed 
that contributions from Glauber gluons cancel in the fi- 
nal cross section, so that we do not need Glauber inter- 
actions in the effective theory. To complete the proof of 
factorization, we now argue that this is indeed the case. 

In principle, Glauber interactions add an additional 
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FIG. 14: RGE running including potential Glauber modes. 



term Cq to the SCET Lagrangian 

£SCET = C na (Xn a ,A a ) + C nb (Xn b ,A s ) + C S {A S ) 

+ C G (A G A.) . (134) 

Glauber interactions in SCET have been considered in 
Refs. [H, H3|, but we will not require an explicit con- 
struction of Cg here. Our arguments will be based on 
the one hand, on the consistency with processes where 
it has been proven that Glauber interactions cancel, and 
on the other hand on systematic scale separation in the 
language of effective field theory. The scale separation is 
valid independently of whether it leads to a factorization 
into simple matrix elements, or whether it leads to a non- 
factorizable matrix element with complicated dynamics. 

The possible danger of the Glauber modes comes from 
the fact that they couple the two collinear sectors n a and 
rib with momentum scaling Q(A 2 ,1,A) and Q(1,A 2 ,A). 
With Cq, there will still be interactions between soft and 
collinear modes present in the Lagrangian even after the 
field redefinition, so we cannot a priori factorize the full 
matrix element into independent soft and collinear ma- 
trix elements. Therefore, we have to revisit each step in 
our derivation with Cg in mind. 

Our argument will be divided into three steps: (i) 
above the scale /is, (ii) at the scale /is, and (hi) below the 
scale fiB- For (i) and (ii) we have to consider Glauber 
modes with momentum scaling Q(A 2 ,A 2 ,A), which we 
call G\ modes. Since they have virtuality ~ /j, B they 
are integrated out at this scale. Any residual effects of 
Glauber interactions below [is could occur from modes 
with momentum scaling Q(A 4 , A 4 , A 2 ), which we call G2. 
These modes arc illustrated in Fig. fTJl 

Above the hard scale [Ih — Q, we have full QCD and 
no distinction between different modes is required, so in 
step (i) we are concerned with contributions of G\ in the 
region fin > fi > \ib- At the scale fin, we integrate out 
hard modes with virtualities Q 2 or higher by matching 
the QCD currents onto SCET currents. For our process, 
the leading operators are given in Eq. (|107p . which con- 
tain only one field per collinear direction. For the theory 



in Eq. (|134p . all other possible operators are power sup- 
pressed. The matching onto these currents is valid at 
an operator level and can be performed with quark and 
gluon states. It is independent of the hadronic matrix 
element we are going to take later on. The key point is 
that the exact same matching calculation and resulting 
Wilson coefficients C occur for threshold Drell-Yan and 
e + e~ — > 2 jets. For these cases it is known 0, HH that 
G\ modes do not affect the matching of the hard func- 
tion H ~ |C| 2 at /in or the running of H in the region 
fiu > fi > fis shown in Figs. [7f b) or^c). The hard 
function H gives a complete description of the physics 
down to the scale \xb whether or not the modes in the 
SCET matrix elements factorize further. In Fig. [TU this 
corresponds to taking the scale /i = [Lb- Therefore, the 
G\ modes can give neither large ln(/xs//ijy) terms nor 
finite contributions above [Lb- 

In step (ii), we integrate out modes with virtualities 
Q 2 A 2 at the scale ps, which may involve matrix elements 
with G\ modes exchanged. This matching affects the n a - 
collinear, nt-collinear, and G\ modes, whose momentum 
scaling below [ib changes to Q(A 4 ,1,A 2 ), Q(1,A 4 ,A 2 ), 
and Q(A 4 , A 4 , A 2 ), respectively. Here we consider the the- 
ory right above fiB including G\ modes, leaving the dis- 
cussion of the theory just below /is and G2 modes to 
step (hi). Thus, we have to consider the matrix element 
of the composite operator 



[Xn a Xn b ]{x + ,x )S(B+ -n a -p a )6(B+ ~n b -p b ) 



[**X»J(0), (135) 



where we suppressed all spin and color indices for sim- 
plicity, and these collinear fields still couple to soft fields 
in their Lagrangians. Since /jb is a perturbative scale, 
we can carry out the matching onto the theory below 
Us at the operator level and do not yet have to consider 
proton states. Since the Glauber gluons are spacelike, 
they cannot cross the final-state cut indicated by the 5 
functions and only appear in virtual subdiagrams. We 
can therefore make a correspondence with the calcula- 
tion in step (i) as follows. For any given final state with 
collinear and soft particles, the SCET computation for 
(ii) is identical to the SCET computation carried out for 
the matching in step (i) but using this particular choice 
of external states. 11 Since that SCET computation can- 
not induce any dependence on G\ in step (i), there can 
also be no contributions from G\ for the forward matrix- 
element computation here. The result of the step (ii) 
matching is thus given by a Wilson coefficient times an 



11 In practice one would never make such a complicated choice, but 
if one does, it must give the same result as picking a minimal 
state for the matching. 
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operator of the form 



J dk+ dk+ C(x + ,x-,B+ - k+,B+ - k+) 

X Xn a (0)T[Yn a Y n b }XnM 6 ( k i ~ *V&) 

x 5{k+ - n b -p b ) 8{u: a - V„ a ) S(uj b - V nb ) 
xx' nb (0)T[YlY na ]x' na (0), 



where the primed collinear fields have scaling Q(A 4 , 1, A 2 ) 
and Q(l, A 4 , A 2 ), and the soft fields in the Y Wilson lines 
have scaling Q(A 2 , A 2 , A 2 ). 

For step (iii) below fis, we have to consider the 
(pp\"'\pp) matrix element of Eq. (|136|) and possible 
contributions from G 2 Glauber gluons, which can now 
also connect to spectator lines in the proton (which are 
primed collinear modes). The G 2 gluons may spoil the 
factorization of the two collinear sectors. To argue that 
this is not the case, we rely heavily on the original proof of 
the cancellation of Glauber gluons for inclusive Drell-Yan 
in Ref. [28|]. By construction, for our observables the k^ b 
variables in Eq. (|136p are of 0(Q\ 2 ) and thus only get 
contributions from the soft gluons. Hence, we are fully in- 
clusive in the Hilbert space of the primed collinear fields. 
Therefore, the G 2 modes as well as possible "ultrasoft" 
Q(A 4 , A 4 , A 4 ) gluons cancel in the sum over states, just as 
in the inclusive case. This discussion for the cancellation 
of G 2 modes is identical to Ref. [29| . where arguments 
were presented for the cancellation of G 2 gluons up to 
the scale induced by the measurement on the final state, 
which in our case is \xb- 

Physically, one could imagine that Glauber modes kick 
the spectators in the proton remnant such that they can 
contribute to B* b . The above arguments show that this 
is not the case, so that our treatment of the proton and its 
remnant in the derivation of the factorization is correct. 

Note that the above arguments do not suffice to show 
that Glauber interactions cancel when there are addi- 
tional hard central jets in the final state. 



D. Final Results for Drell-Yan 

In this subsection, we present the final results for the 
isolated Drell-Yan cross section. Our discussion is split 
into four parts: the leptonic tensor, the hard function, 
the soft function, and the final cross section for beam 
thrust. 



1. The Leptonic Tensor 

To give an explicit example, we now apply the final 
factorization result in Eq. (|132[) to the Drell-Yan process 
with L = £ + £~. The relevant QCD currents are the 
vector and axial-vector currents Jj^j, with h — {V, A}, 
already given in Eq. (|98p . The corresponding leptonic 



contributions arc 

47ra e m 



r 



-Qfu(p 2 h fJ, v(j> 1 ) 



+ 



1 - m z /q z 



(136) L l X f {p u p 2 ) 



47ra r 



q 2 1 — m 2 z I q 2 
X w(p 2 )7 M (^ - aab)v{Pi) , 



(137) 



where in this subsection p\ = pi+ and p 2 = pe- are the 
lepton momenta, Qf is the quark charge (in units of |e|), 
and vij, aij are the standard vector and axial couplings 
of the leptons and quarks of flavor / to the Z boson. We 
will include the width of the Z later in Eq. (|147jl . 

The leptonic phase space integral is 



/ 



d 4 p 1 d 4 p 2 xf 2 , x( 2u4 / _n a + n b 
—j—^—5( Pl )5(p 2 )6 [q —+q + —-pi-p 2 



1 



dAy dip 
32 7T 2 J 1 + cosh Ay 



(138) 



where ip is the azimuthal angle of the leptons in the trans- 
verse plane and Ay is the rapidity difference of the two 
leptons: 



j/i = i In H±Pl 5 Ay = y l -y 2 . (139) 
2 n a -pi 



Since we expanded <jV = 0, the leptons are back-to-back 
in the transverse plane, which implies that at the order 
we are working 



— = — = —, Y = - (?/i + y 2 ) ■ 
Pi Pi 1 2 



Pit 



-P 2 t , 



-2 -e 

Pit — Pit 



2(1 + cosh Ay) 



(140) 



Thus, the leptonic kinematics is described by the four 
independent variables {q 2 , Y, Ay, ip}, with {Y, Ay} being 
equivalent to {2/1,2/2}- For simplicity, we assume that 
we do not distinguish the two leptons, as one would for 
example by measuring their rapidities yi or transverse 
momenta pit- We can then integrate over < tp < 2-k 
and —00 < Ay < 00 in Eq. (|138() . giving an overall factor 
of 47r. The leptonic tensor, Eq. (|101|) . now becomes 



1 



dAy dip 
32tt 2 J 1 + cosh Ay 



L hf(Pi>P2) L h>f'(PiiP2) 



spins 



iq 2 



(±*- g ^) Lhh , tf ,{ q *). 



(141) 
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where 



L vvff ,{q 2 ) = Q f Qf 



(QfVf> +v f Qf,)vt 



+ 



1 — m z 1 1 q 2 
v f v r {vj+a 2 ) 



LAAff(q 2 



L A vfr(q 2 ) 



CLfClf^v] + of) 

(l-m|/g 2 ) 2 
~ a f 



1 - m 2 z /q 2 
L VA f>f{q 2 ) 



Qfve 



v r (v 2 + a 2 e )~l 



1 - m 2 z /q 2 



(142) 



2. The Hard Function 



Using parity and charge conjugation invariancc of 
QCD, the matching coefficients for the vector and axial- 
vector QCD currents can be written as 

Cyf qq (b a ,b b ) = C vfq (q 2 )( 1 ^, 

C^(baM) = C Afq (q 2 )(^ ±l5 re, 

C^ikM) = C Ag (q 2 ) (b a + b b yie<>° XK blb K b . (143) 



By Lorentz invariance (or reparametrization invariancc 
of n a _ b and n a ^ b [86|). the scalar coefficients can only de- 
pend on b a ■ b b = x a x b E 2 m = q 2 . In principle, parity 
and charge conjugation would also allow the Dirac struc- 
tures (b a - hy 5 af3 and (b a - b b y (7 5 ) Q/3 . However, as 
the vector and axial- vector currents are chiral even and 
the matching from QCD conserves chirality for massless 
quarks, these cannot be generated. For the gluon opera- 
tor, the symmetry of the Wilson coefficient [see Eq. (11081) ] 
requires it to be proportional to q^ — b^ + b^. Current 
conservation for the vector current requires q^Cy j qq - = 0, 
which eliminates this term. Thus, as expected, the only 
contribution for the gluon operator is due to the axial 
anomaly, coming from the diagram in Fig. |6(b)| Since 
we neglect the lepton masses, q^L^j = 0, and thus 
C A f g g does not survive the contraction of the leptonic 
and hadronic tensors for L = . Hence, the gluon 
beam functions do not contribute to Drell-Yan, and the 
gluon PDF only appears through its contribution to the 
quark beam functions. Inserting Eq. (|143jl into the gen- 
eral expression for the hard function in Eq. (|125p , we 
obtain 



H hh'ff q <j(ba,b b ) 



H 



hh'ff'qq&aM) 



1 

1 . 
if 

4 AT, 



<T - \{« + «)] C* hfq (q 2 )C h , f/q (q 2 ) 
^x K n x a nSC* hfq (q 2 )C h , f , q (q 2 ) 



(for hti = {VV, AA}) , 
(for hti = {VA,AV}). 



(144) 



At one loop, the vector and axial-vector coefficients are equal and diagonal in flavor and the SCET matching 
computation was performed in Refs. (8tj l87j. in agreement with the one- loop form factors 



Cv Sq {q 2 ) = C Afq (q 2 ) = S fq C(q 2 ) , C{q\ fi) = 1 + 



47T 



In' 



it) 



3 In 



it) 



7T 

~6~ 



(145) 



The vector current coefficient at two loops was obtained in Refs. [121 |88| from the known two-loop quark form 
factor (89W92I ]. Starting at three loops, it can have a contribution that is not diagonal in flavor, i.e., is not proportional 
to 5f q . The axial- vector coefficient can also receive additional diagonal and nondiagonal contributions starting at two 
loops from the axial anomaly [93l - [95| . The anomaly contributions cancel in the final result in the sum over / as long 
as one sums over massless quark doublets. Therefore, they will cancel when the hard matching scale is much larger 
than the top-quark mass, in which case the top quark can be treated as massless. On the other hand, they have to 
be taken into account when the matching scale is below the top-quark mass, in which case the top quark is integrated 
out during the matching step and its mass cannot be neglected. 

Combining Eqs. ([Til)) and (flUj)) with Eqs. (flUj) and ([US]) , the coefficients H i:j (q 2 ,Y) in Eq. ([T3"3")) become 



1 



2E 2 



3q 2 



-E 

N °fp 



iw/f(9 2 )C^(g 2 )Cv f ,( 9 2 ) + L^/Hg 2 )^,!? 2 )^,^)] ^ <r ff M -( g 2 ^) , (146) 



where at one loop 



47TOC 



CO 



3N c E 2 cm q 2 ' 



Q 



(v 2 + a 2 q )(v 2 + of) - 2Q q v q v e (l - m 2 z /q 2 ) 
(l-m|/ 9 2)2 + m |r|/g4 



\C(q\n)\\ (147) 
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with \C(q 2 , fi)\ 2 given by Eq. (j 145[) . and where we also included the nonzero width of the Z . The RGE for the hard 
function H qq (q 2 ,fi) is 



— =lH(q 2 ,v)Hqq(q 2 ,^), 7#(<? 2 ,m) = 2r cusp [a s (/i)]ln^2 +7j?[a s (^)] , 



d/i * fi 

where T cusp is the universal cusp anomalous dimension [74j], and the one-loop non-cusp term is 7#[a s (//) 
-3a s (/i) C F /ir H3- The solution of Eq. (fT48|) has the standard form 

- 2 \ Vh(ho,iA 



(148) 



H gg -(q 2 ,fl)=H qg -(q 2 ,fl )U H (q 2 > fl (h fl), U H tf, fl , fi) = e K "^ (%) 



where Knifio^fi) and t\h (a*o,m) are analogous to Eq. ([63 



KH{flQ,fl) 



da s 



4 r CUS p(o; s ) 



da' 



i (Mo) 



/3(a 



»W da. 



i (Mo) 



(149) 



r cusp (a s ). (150) 



Together, Eqs. (|149| and ()150|) sum the large logarithms occurring in isolated Drcll-Yan between the scales ixh and 
fiB- Elcctroweak corrections to the hard function H qq (q 2 , fi) can be included using the results of Refs. 



3. The qq Soft Function 

The incoming hemisphere soft function contains incoming Wilson lines stretching from — oo to along n a and rib- 
Under time reversal, each incoming Wilson line transforms into a corresponding outgoing Wilson line stretching from 
to oo along the opposite direction, 



T _1 K T = Pexp 



ig / dsn b -A s (sn b ) 



(i 



Yn b , (151) 



where P denotes anti-path ordering. Since T itself does not affect the original ordering of the field operators, time 
ordering turns into anti-time ordering and vice versa. In addition Tn a f^ s T _1 = n b -pb iS . Therefore, time-reversal 
invariance implies 

S$Li(K, K) = ^ tr {0\T[YlY na ]S(k+ - n b -p b , s ) S(k+ - n a -p a , s )T[Y r \Y nb ] |0)* 

= i-tr (0\T[YlY nb ]6(k+ - n a -p a>s ) S(k+ - n b -p b , s )T[Y^ b Y na ] |0) . (152) 

In the second step, the complex conjugation has no effect since the matrix element is real, and we used parity 
to switch n bt a back to n a%b . For comparison, the hemisphere soft function with o utgoing Wilson appearing in the 
double-differential hemisphere invariant-mass distribution in e + e~ — > 2 jets [69l. l70l. [98l4lOl| is 

S$Li{k+, k+) = i- tr (0\T[YlY nb ] 5(k+ - n a -p a , s ) 5(k+ - n b -p b , s ) T^Y^] |0) . (153) 

This is almost the same as Eq. (|152[) . the only difference being the opposite time ordering. Thus, Sihcmi and Shemi 
are equal at one loop, where the time ordering is still irrelevant. Beyond one loop, S'ihemi and S^emi may in general be 
different. However, since the beam and jet functions have the same anomalous dimension, the combined anomalous 
dimension of the hard and beam functions in isolated Drcll-Yan agrees with that of the hard and jet functions for 
the e + e~ hemisphere invariant- mass distribution. The consistency of the RGE in both cases then requires that Sihcmi 
and Shemi have the same anomalous dimension to all orders in perturbation theory. In addition, the purely virtual 
contributions, obtained by inserting the vacuum state, are the same in both cases, 

«f ual (fc Q + , K) = ^ KK) 5(ki)tr | (0\T[YlY nb ] | ) | 2 = S^t^, *+) . (154) 

Hence, S'ihcmi and Shemi can only differ by finite real- emis sion corrections at each order in perturbation theory. 
Using the one- loop results for S^emi f rom Refs. (70l llOC| . we have 

S!l mi (kt, K) = 6{kt) 5{kt) + S(k+) 5 llo °P(fc+) + 5 llo °P(fc+) 5(k+) , 

~B{k+/fi) ln(fc+//i) 



Am \ fi 



k+/fL 



Y d(k + )}. (155) 



The plus distribution is defined in Eq. (|59p . The one-loop soft function for beam thrust in Eq. (|2l)|) then becomes 
S B (k + ,fi) =S(k+) + 2S llo °v(k + ). 
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4- Final Cross Section for Beam Thrust 
The differential cross section for beam thrust in Eq. (|32|) including the RGE running is 



da 



dq 2 c\YAt b 



a Q 'Y^H lJ (q 2 ,fi H )UH(q 2 ,fJ-H,tJ-s) J dt a dt b Q S b (Qt b - 



Q 



x / dt' a Bi(t a — t' a , x a , Hb) U B (t' a , fJ>B,fJ>s) j dt b Bj(t b - t' b ,x b , n B )U B (t b , lib, Us) ■ (156) 

For simplicity, we evolve the hard and beam functions from their respective hard and beam scales, \iu and fi Bl down 
to the common scale fi = fis of the soft function. In this way, we do not need to consider the running of the soft 
function separately. Different choices for /i are all equivalent, as we discussed in Sec. Ill Ft At LL, we include the 
one- loop cusp anomalous dimension in the evolution kernels U B and U B , and at NLL we include the two-loop cusp 
and one-loop non-cusp anomalous dimensions. In both cases we use the LO results as initial conditions. 

We also consider the fixed-order a s expansion. To our knowledge dcr/dg 2 dydrB has not been considered in per- 
turbation theory in full QCD even at one loop. To obtain an expression for da/dq 2 dYdr B at NLO in a s and leading 
order in the power counting, we drop the evolution factors U B and U B and expand all functions to NLO at a common 
scale 11. From the above NLO results for the hard and soft functions and the NLO results for the beam functions 
from Sec. IIII1 we find 



da 



dq 2 dYdr B 



E 



Q 



2 , K 2 + of ){vj + of) - 2Q l v i v e (l - m 2 z /g 2 ) 
(l-m%/q*) 2 + m%r%/q4 



f d£a d& [X a X b 2 \ 



(157) 



Here, /i/ a (£a,A*) and fj/b(£,b, fJ-) are the PDFs for parton i in proton a and parton j in (anti-)proton b. At tree level, 
the nonzero coefficients are 



C*™ e (0 a ,z 6 , g 2 ,T B ,/i) = C q r q cc {z a ,z b ,q 2 ,T B ,^) = S(t b )S(1 - z a ) 5(1 - z b ) . 
At one loop, we obtain 



(158) 
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(159) 



The coefficient C gg only starts to contribute at two loops. 
The single logarithms of q 2 / /j, 2 are multiplied by the QCD 
splitting kernels and are resummed by the PDFs. Thus, 
in fixed-order perturbation theory the PDFs should be 
evaluated at the hard scale fj, = Q, such that there are 



no large logarithms when integrating over < t b < 1. 
However, if the integration is restricted to t b < r B ut <C 
1, the plus distributions in t b produce large logarithms 
In 2 t™' and ln.T B ut , which make a fixed-order expansion 
unreliable. These are precisely the logarithms that are 
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rcsummcd by the combined RGE of hard, jet, and soft 
functions in Eq. (|156jl . 

V. ISOLATED DRELL-YAN CROSS SECTION 

In this section we illustrate our results with plots of the 
isolated Drell-Yan cross section. Rather than considering 
the cross section as a function of two variables, B+ and 
Bt , we use the beam thrust tb defined in Sec. Ill Bl We 
consider both the differential cross section da/dQclYdrs 
as a function of tb, see Eq. (j3"2")h as well as the cross 
section integrated over < tb < cxp(— 2y™*) as a func- 
tion of y c B ut , see Eqs. (EU) and (J33]). We fix Q = y/q 1 
to a few representative values. We also restrict our dis- 
cussion to back-to-back leptons, Y — 0, since measuring 
tb mainly affects the normalization and not the shape of 
the rapidity distribution. 

For our cross-section predictions, we always use the 
NLO PDFs from MSTW2008 [n!, with the correspond- 
ing a s (mz) = 0.1202 and two- loop five-flavor running for 
a s (fi). 

We show results for the cross section both in the fixed- 
order expansion at LO and NLO [sec Eq. f| L57[) ] and in the 
rcsummcd expansions at LL and NLL [see Eq. (|156p ]. For 
the resummed results we always choose the hard, beam, 
and soft scales as — fi, /is — H^/tb, fJ-s = ^ t b, 
and for the fixed-order results we use a common fixed 
scale /i. The central values correspond to fx = Q, and 
the bands show the scale variation obtained by varying 
/.t = 2Q and /i = Q/2. Since our purpose here is to 
illustrate the main features of the factorized cross section 
with beam functions, we will limit ourselves to the LL and 
NLL results without also including additional fixed-order 
corrections to the hard, beam, and soft functions. A 
complete analysis combining both NLO corrections and 
NNLL resummation, and a more detailed analysis of scale 
uncertainties, is left for future work. 

In Fig.QjjJ we show results for pp collisions at the LHC 
with E cm = 7 TcV and two different values Q = 100 GeV 
and Q = 1 TeV. For Y = this corresponds to x a = xj, = 
0.014 and x a = Xb = 0.14 respectively, representing two 
typical values for the isolated factorization theorem. In 
Fig. [TH we show results for pp collisions at the Tevatron 
for Q = 500 GeV corresponding to x a = xt, = 0.26. 

The top row in Fig.Ql)]depicts the LHC cross section as 
a function of tb at several different orders. At tree level, 
there is only a 5 function and no radiation. The green 
(light) curve and band show the NLO-expanded cross 
section, which grows for decreasing tb , showing that the 
radiation is peaked in the forward direction. As tb — >• 
it exhibits the expected singular behavior, and this IR 
singularity is canceled in the integrated cross section by 
a corresponding S function at tb = 0. The blue (dark) 
and orange (medium) curves are the resummed results at 
LL and NLL, respectively. As expected, the resumma- 
tion has a large effect at small tb and effectively regulates 
the IR singularity in the NLO result. The curves are not 



plotted for t b < 0.007 in the left panel and t b < 0.001 
in the right panel, because at this point the soft scale 
drops below 1 GeV. Near this cutoff, the soft function 
becomes nonperturbative and so our purely perturbative 
results should not be taken too seriously. Nevertheless, 
it is interesting to sec that the resummed results show a 
characteristic turnover that one would also expect from 
nonperturbative corrections. These nonperturbative cor- 
rections scale as Aqcd /(Qtb) and hence become less rel- 
evant than the perturbative corrections for tb > 0.02 
and tb > 0.002 in the left and right panels, respectively. 

In the middle row of Fig. [15j wc compare our results 
expanded to fixed-order at LO (gray horizontal band) 
and NLO (green band with solid central line) as a func- 
tion of y B ut - Since we plot the integrated cross section, 
the tree-level S function gives a constant contribution in 
y B ut - As expected, the NLO corrections to the LO re- 
sult become very large as we move away from y B nt = 
because of the large double logarithms in the fixed-order 
expressions. Hence, a fixed-order expansion is not reli- 
able here. The green dashed line illustrates what hap- 
pens if we take the NLO result given by the solid line 
but exclude the contribution of the gluon PDF to the 
quark beam function. The gluon contribution has a big- 
ger effect at smaller Q, because it corresponds to smaller 
x. Overall, it reduces the integrated cross section, which 
is the same effect we already observed in the integrated 
beam function in Sec. Mil The plots in the bottom row 
of Fig. [15] show the LL, NLL, and NLO results, which 
are the integrated versions of the curves in the top row. 
Here, larger y c B ut towards the left of the plot corresponds 
to smaller tb . A large y B ut implies a stronger constraint 
on the final state and hence a smaller cross section. In 
this region, the logarithmic resummation is necessary and 
suppresses the cross section. For y c B ut —> the LL and 
NLL approach the LO result as they must. 

Figure [1^] shows the corresponding plots for the Teva- 
tron. The left panel and right panel are equivalent to the 
top row and bottom row of Fig. [T5j respectively showing 
the differential and integrated cross sections. The plots 
show similar features overall, though the overall scale un- 
certainties are slightly larger here, because the PDFs are 
evaluated at a larger x. In the Tevatron case, even though 
we are at larger x compared to the right panels in Fig. 1151 
the cross section is larger due to the presence of valence 
antiquarks (note the different units in the two plots). 



VI. CONCLUSIONS 

Experimental measurements at the LHC or Tevatron 
are typically characterized by two conditions. First, the 
dominant part of the cross section arises from parton mo- 
mentum fractions x away from one, and second, to probe 
the hard scattering the measurements impose restrictions 
on the final state to identify and isolate hard leptons or 
jets. 

Factorization is required to separate the perturbatively 
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FIG. 15: Cross sections for beam thrust at the LHC with E cm = 7 TeV at Y = and Q = 100 GeV (left column) and Q = 1 TeV 
(right column). Top row: The cross section differential in tb at NLO, LL, and NLL. Middle and bottom rows: The cross section 
integrated up to tb < exp(-2j/g ut ) at LO, NLO, LL, and NLL. 



calculable pieces from the nonperturbative parton distri- 
bution functions, and is a key ingredient for the resum- 
mation of large logarithms that occur due to phase-space 
restrictions. The most well-known factorization theo- 
rem for inclusive Drell-Yan applies for generic momen- 
tum fractions, but the hadronic final state is completely 
summed over, only subject to overall momentum con- 
servation. This requires an inclusive experimental mea- 



surement, with only mild restrictions on the final state. 
On the other hand, threshold factorization theorems for 
Drell-Yan or dijet production take into account phase- 
space restrictions and resum resulting large logarithms, 
but they are only valid in the limit x — > 1. Thus each of 
these cases satisfies only one of the above experimental 
conditions. 

In this paper we have studied factorization for generic x 
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and with explicit restrictions on the hadronic final state. 
We considered the simplest situation, namely Drell-Yan 
production pp — > X£ + £~ at generic x with a restriction 
on X that vetoes hard central jets, which we call isolated 
Drell-Yan. The restriction on the hadronic final state is 
implemented by dividing the total hadronic momentum 
into two hemispheres and requiring the components _B+ 
and B b of the resulting hemisphere momenta to be small. 
For this situation, we prove a factorization theorem for 
the cross section differential in the hadronic variables 
and B b (as well as the leptonic phase-space variables) . It 
allows us to systematically resum large double logarithms 
of B^ b /Q arising from the phase-space restrictions. The 
factorization theorem also applies to other isolated pro- 
cesses of the form pp — > XL, like H 77 or H —> A£. 

The main conclusion from our analysis is that PDFs 
alone are insufficient to properly describe the initial state 
in the collision. The restriction on the hadronic final state 
effectively probes the proton prior to the hard interaction 
by constraining the virtuality of the colliding hard par- 
ton. At that point, the parton cannot be confined to 
the proton anymore and must properly be described as 
part of an incoming initial-state jet. This description is 
given by quark and gluon beam functions, which replace 
the PDFs in the factorization theorem. In addition to 
the parton's momentum fraction, the beam functions de- 
pend on the virtuality of the colliding parton and live at 
an intermediate beam scale /is of the order of this vir- 
tuality. At the scale fi = /is, the beam functions can 
be matched onto PDFs evaluated at this /1 times per- 
turbatively calculable corrections. In this way, our fac- 
torization theorem unambiguously determines the proper 
scale at which the PDFs must be evaluated. This scale 
Us is much smaller than the hard scale of the partonic 
collision. As a result, the evolution of the initial state 
below [is is governed by the PDFs, while above /is it 
is governed by the beam function, whose evolution de- 
pends on the parton's virtuality rather than momentum 



fraction. Other differences compared to the PDF evolu- 
tion are that the beam-function evolution sums double 
logarithms and does not involve parton mixing. 

In a Monte Carlo setting the corresponding physical 
effects should be described by the initial-state parton 
shower. Our factorization theorem thus provides a way to 
explicitly check whether the initial-state parton shower 
rcsums the correct double logarithms, which is left for 
future work. We believe that experimental measure- 
ments of the isolated Drell-Yan spectrum will provide 
a direct method for testing the initial-state shower in 
Monte Carlo. This spectrum combined with the results 
reported here, is therefore useful for tuning the Monte 
Carlo with early LHC data. 

Even though our derivation is only rigorous for the 
specific case of isolated Drell-Yan, we argued that the 
necessity for beam functions is more general, essentially 
applying to any process where the hadronic final state is 
restricted in a similar way. For example, the consistency 
of the renormalization group evolution implies that the 
description for any threshold process can be extended to 
a respective isolated case by supplementing it with cor- 
responding variables B^ b , replacing the PDFs by beam 
functions, and replacing the threshold soft function by 
an appropriate isolated soft function. 

We briefly discussed the extension of isolated Drell-Yan 
to isolated dijet production. We also pointed out that if 
the hadronic final state is constrained by different global 
variables, one can expect to find different beam functions 
to encode these constraints, which will then sum double 
logarithms in these variables. 

Ultimately, we hope that this type of factorization the- 
orem with beam functions will bridge the gap between 
experimentally realistic cuts for LHC measurements and 
systematically improvable theoretical results that go be- 
yond fixed-order calculations. 
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